A Floating Point Package 
for the HPC 



INTRODUCTION 

This report describes the implementation of a Single Preci- 
sion Floating Point Arithmetic package for the National 
Semiconductor HPC microcontroller. The package is based 
upon the IEEE Standard for Binary Floating-Point Arithmetic 
(ANSI/IEEE Std 754-1985). However, the package is not a 
conforming implementation of the standard. The differences 
between the HPC implementation and the standard are de- 
scribed later in this report. 

The following single precision (SP) operations have been 
implemented in the package. 

(1) FADD. Addition of two SP floating point (FLP) numbers. 

(2) FSUB. Subtraction of two SP FLP numbers. 

(3) FMULT. Multiplication of two SP FLP numbers. 

(4) FDIV. Division of two SP FLP numbers. 

(5) ATOF. Convert an ASCII string representing a decimal 

FLP number to a binary SP FLP number. 

(6) FTOA. Convert a binary SP FLP number to a decimal 

FLP number and output the decimal FLP number 

as an ASCII string. 
The report is organized as follows. The next section dis- 
cusses the representation of FLP numbers. Then, the differ- 
ences between the HPC implementation and the IEEE/ 
ANSI standard are described. This is followed by a descrip- 
tion of the algorithms used in the computations. Appendix A 
is a User's Manual for the package, Appendix B describes 
the test data for the package and Appendix C is a listing of 
the code. 

Note that this report assumes that the reader is familiar with 
the IEEE/ANSI Binary Floating-Point Standard. Please refer 
to this document for an explanation of the terms used here. 

REPRESENTATION OF FLOATING POINT NUMBERS 

The specification of a binary floating point number involves 
two parts: a mantissa and an exponent. The mantissa is a 
signed fixed point number and the exponent is a signed 
integer. The IEEE/ANSI standard specifies that a SP FLP 
number shall be represented in 32 bits as shown in Figure 1 . 
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FIGURE 1 
The significance of each of these fields is as follows: 

1.S — this 1-bit field is the sign of the mantissa. S = 
means that the number is positive, while S = 1 means 
that it is negative. 

2. E — this is the 8-bit exponent field. The exponent is repre- 
sented as a biased value with a bias of 127-decimal. 

3. F — this is the 23-bit mantissa field. For normalized FLP 
numbers (see below), a MSB of 1 is assumed and not 
represented. Thus, for normalized numbers, the value of 
the mantissa is 1.F. This provides an effective precision 
of 24 bits for the mantissa. 

Normalized FLP number: A binary FLP number is said to be 
normalized if the value of the MSB of the mantissa is 1. 
Normalization is important and useful because it provides 
maximum precision in the representation of the number. If 
we deal with normalized numbers only (as the HPC imple- 
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mentation does) then since the MSB of the mantissa is al- 
ways 1, it need not be explicitly represented. This is as 
specified in the IEEE/ANSI standard. 
Given the values of S , E and F, the value of the SP FLP 
number is obtained as follows. 

If < E < 255, then the FLP number is (-1) 

»S*1.F*2*(E-127). 

If E = 0, then the value of the FLP number is 0. 

If E = 255, then the FLP number is not a valid number 

(NAN). 
The above format for binary SP FLP numbers provides for 
the representation of numbers in the range -3.4*10*38 to 
-1.75*10»-38, 0, and 1.75*10*-38 to 3.4*10*38. The ac- 
curacy is between 7 and 8 decimal digits. 

DIFFERENCES BETWEEN THE IMPLEMENTATION AND 
THE IEEE/ANSI STANDARD 

The IEEE/ANSI standard specifies a comprehensive list of 
operations and representations for FLP numbers. Since an 
implementation that fully conforms to this standard would 
lead to an excessive amount of overhead, a number of the 
features in the standard were dropped. This section de- 
scribes the differences between the implemented package 
and the standard. 

1. Omission of -0. The IEEE/ANSI standard requires that 
both + and — zero be represented, and arithmetic car- 
ried out using both. The implementation does not repre- 
sent — 0. Only +0 is represented and arithmetic is car- 
ried out with + only. 

2. Omission of Infinity Arithmetic. The IEEE/ANSI standard 
provides for the representation of plus and minus Infinity, 
and requires that valid arithmetic operations be carried 
out on Infinity. The HPC implementation does not support 
this. 

3. Omission of Quiet NaN. The IEEE/ANSI standard pro- 
vides for both quiet and signalling NaNs. The HPC imple- 
mentation provides for signalling NaNs only. A signalling 
NaN can be produced as the result of overflow during an 
arithmetic operation. If the NaN is passed as input to fur- 
ther floating point routines, then these routines will pro- 
duce another NaN as output. The routines will also set 
the Invalid Operation flag, and call the user floating point 
error trap routine at address FPTRAP. 

4. Omission of denormalized numbers. Denormalized num- 
bers are FLP numbers with a biased exponent, E of zero 
and a non zero mantissa F. Such denormalized numbers 
are useful in providing gradual underflow to zero. Denor- 
malized numbers are not represented or used in the HPC 
implementation. Instead, if the result of a computation 
cannot be represented as a normalized number within the 
allowable exponent range, then an underflow is signaled, 
the result is set to zero, and the user floating point error 
trap routine at address FPTRAP is called. 

5. Omission of the Inexact Result exception. The IEEE/ 
ANSI standard requires that an Inexact Result exception 
be signaled when the rounded result of an operation is 
not exact, or it overflows without an overflow trap. This 
feature is not provided in the HPC implementation. 
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6. Biased Rounding to Nearest. The IEEE/ANSI standard 
requires that rounding to nearest be provided as the de- 
fault rounding mode. Further, the rounding is required to 
be unbiased. The HPC implementation provides biased 
rounding to nearest only. An example will help clarify this. 
Suppose the result of an operation is .b1b2b3XXX and 
needs to be rounded to 3 binary digits. Then if XXX is 
OYY, the round to nearest result is .b1b2b3. If XXX is 
1 YY, with at least one of the Y's being 1 , then the result is 
.b1 b2b3 + 0.001 . Finally if XXX is 1 00, it is a tie situation. 
In such a case, the IEEE/ANSI standard requires that the 
rounded result be such that its LSB is 0. The HPC imple- 
mentation, on the other hand, will round the result in such 
a case to .b1b2b3 + 0.001. 

DESCRIPTION OF ALGORITHMS 

1 . General Considerations. The HPC implementation of 
the SP floating point package consists of a series of sub- 
routines. The subroutines have been designed to be com- 
patible with the CCHPC C Cross Compiler. They have, 
however, not been tested with the CCHPC Cross Compil- 
er. 

The Arithmetic subroutines that compute F1 op F2 
(where op is +, -, * or /) expect that F1 and F2 are 
input in the IEEE format. Each of F1 and F2 consists of 
two 16-bit words organized as follows. 
Fn-HI: S EXP 7 MS bits of F 
Fn-LO: 16 LS bits of F 

In the above, S is the sign of the mantissa, EXP is the 
biased exponent, and F is the mantissa. 
On input it is assumed that F1 -HI is in register K, F1 -LO 
is in the accumulator A, and F2-HI and F2-LOareon the 
stack just below the return address i.e., F2-HI is at 
W(SP-4) and F2-LO is at W(SP-6). The result, C, is also 
returned in IEEE format with C-HI in register K and C-LO 
in the accumulator A. 

The two Format Conversion routines, ATOF and FTOA 
expect that on entry, register B contains the address of 
the start of the ASCII byte string representing the decimal 
FLP number. ATOF reads the byte string starting from 
this address. Note that the string must be terminated with 
a null byte. The binary floating point number is returned in 
registers K and A. FTOA, on the other hand, writes the 
decimal FLP string starting from the address in register B 
on entry. A terminating null byte is also output. Also, 
FTOA expects that the binary FLP number to be convert- 
ed is in registers K and A on entry. 
Most of the storage required by the subroutines is ob- 
tained from the stack. Two additional words of storage in 
the base page are also used. The first is W(0), and is 
referenced in the subroutines as W(TMP1). The second 
word of storage can be anywhere in the base page and is 
used to store the sticky flags used to signal floating point 
exceptions. This is referenced in the subroutines as 
W(FPERWD). Thus any user program that uses the float- 
ing point package needs to have the symbols TMP1 and 
FPERWD defined appropriately. 

2. Exception Handling. The following types of exception 
can occur during the course of a computation. 

(i) Invalid Operand. This exception occurs if one of the 
input operands is a NaN. 

(ii) Exponent Overflow. This occurs if the result of a com- 
putation is such that its exponent has a biased value 
of 255 or more. 



(iii) Exponent Underflow. This occurs if the result of a 
computation is such that its exponent is or less. 

(iv) Divide-by-zero. This exception occurs if the FDIV rou- 
tine is called with F2 being zero. 

The package signals exceptions in two ways. First a word 
at address FPERWD is maintained that records the histo- 
ry of these exception conditions. Bits 0-3 of this word are 
used for this purpose. 
Bit — Set on Exponent Overflow. 
Bit 1 — Set on Exponent Underflow. 
Bit 2— Set on Illegal Operand. 
Bit 3 — Set on Divide-by-zero. 

These bits are never cleared by the floating point pack- 
age, and can be examined by the user software to deter- 
mine the exception conditions that occurred during the 
course of a computation. It is the responsibility of the user 
software to initialize this word before calling any of the 
floating point routines. 

The second method that the package uses to signal ex- 
ceptions is to call a user floating point exception handler 
subroutine whenever an exception occurs. The corre- 
sponding exception bit in FPERWD is set before calling 
the handler. The starting address of the handler should 
be defined by the symbol FPTRAP. 

3. Unpacked Floating Point Format. The IEEE/ANSI stan- 
dard floating point format described earlier is very cum- 
bersome to deal with during computation. This is primarily 
because of the splitting of the mantissa between the two 
words. The subroutines in the package unpack the input 
FLP numbers into an internal representation, do the com- 
putations using this representation, and finally pack the 
result into the IEEE format before return to the calling 
program. The unpacking is done by the subroutine FUN- 
PAK and the packing by the subroutine FPAK. The un- 
packed format consists of 3 words and is organized as 
follows. 

Fn-EXP.Fn-SIGN 8 bits biased sign (extended to 
exponent 8 bits) 

Fn-HI MS 16 bits of mantissa 

(implicit 1 is present as MSB) 

Fn-LO LS 8 bits of Eight 

mantissa Zeros 

Since all computations are carried out in this format, note 
that the result is actually known to 32 bits. This 32-bit 
mantissa is rounded to 24 bits before being packed to the 
IEEE format. 

4. Algorithm Description. All the arithmetic algorithms first 
check for the easy cases when either F1 or F2 is zero or 
a NaN. The result in these cases is immediately available. 
The description of the algorithms below is for those cases 
when neither F1 nor F2 is zero or a NaN. Also, in order to 
keep the algorithm description simple, the check for un- 
derflow/overflow at the various stages is not shown. The 
documentation in the program, the descriptions given be- 
low, and the theory as described in the references should 
allow these programs to be easily maintained. 

(i) FADD. 

The processing steps are as follows: 
1. Compare F1-EXP and F2-EXP. Let the difference 
be D. Shift right the mantissa (Fn-HI. Fn-LO, n = 1 
or 2) of the FLP number with the smaller exponent 
D times. Let the numbers after this step be F1- 
EXP.F1-SIGN, F1-HI, F1-LOand F2-EXP.F2-SIGN, 



F2-HI and F2-LO. This step equalizes the two expo- 
nents. 

2. Take the XOR of F1 -SIGN and F2-SIGN. If this is 0, 
then go to step 4, else go to step 3. 

3. Do a true subtract of F2-LO from F1-LO. (A true 
subtract is when the SUBC instruction is preceded 
by a SET C instruction.) Then do a 1 's complement 
subtract of F2-HI from F1-HI. If the last subtract 
resulted in C = 1, then go to step 3.2, else go to 
step 3.1. 

3.1. Get here means that F2 is larger than F1, and 
the computed result is negative. Take the 2's com- 
plement of the result to make it positive. Set the 
sign of the result to be the sign of F2. Go to step 
3.3. 

3.2. Get here means F1 is larger than F2, and the 
result of the mantissa subtract is positive. Set the 
sign of the result to be the sign of F1 . Go to step 
3.3. 

3.3. The result after a subtract need not be normal- 
ized. Shift left the result mantissa until its MSB is 1 . 
Decrement the exponent of the result by 1 for each 
such left shift. Go to step 5. 

4. Add F2-LO to F1-LO. Next add with any carry from 
the previous add, F2-HI to F1-HI. If this last add 
results in C = 1, then go to step 4.1, else go to 
step 5. 

4.1 . Rotate Right with carry C-HI. Next load C-LO in 
and rotate it right with carry. Increase the exponent 
of the result, C by 1 . Go to step 5. 

5. Round the result. Go to step 6. 

6. Pack the result and return, 
(ii) FSUB. 

The processing steps are as follows: 

1 . Copy F2 to the stack and change its sign. Go to 
step 2. 

2. Call FADD. 

3. Remove the copy of -F2 from the stack and return, 
(iii) FMULT. 

The processing steps are as follows. 

1. Add F1-EXP and F2-EXP to get C1-EXP. Subtract 
from C1-EXP 127-decimal which is the IEEE bias, 
to get C-EXP. Go to step 2. 

2. Take the XOR of F1-SIGN and F2-SIGN to get C- 
SIGN. Go to step 3. 

3. Compute F1-HI*F2-HI. Let the upper half of the 
product be C1-HI and the lower half C1-LO. Go to 
step 4. 

4. Compute F1-HI*F2-LO. Let the upper half of this 
product be C2-HI. Add C2-HI to C1-LO to give 
C11-LO. If this last add results in C = 1, then 
increment C1-HI. Go to step 5. 

5. Compute F1-LO*F2-HI. Let the upper half of this 
product be C3-HI. Add C3-HI to C11-LO to get 
C12-LO. If this last add results in C = 1, then 
increment C1-HI. Go to step 6. 

6. Mantissa normalization. If the MSB of C1-HI is 1, 
then increment C-EXP, else shift left C1-HI.C12- 
LO. Go to step 7. 

7. Round C1-HI.C12-LO to get C-HI.C-LO. Go to 
step 8. 



8. Pack C-EXP.C-SIGN, C-HI and C-LO and return as 
the answer, 
(iv) FDIV. 

The processing steps are as follows: 

1. Compare F1 -HI and F2-HI. If F2-HI is greater than 
F1-HI then go to Step 3, else go to step 2. 

2. Shift right F1-HI.F1-LO. Increase F1-EXP by 1. 

3. Subtract F2-EXP from F1-EXP. Add to the result 
127-decimal to get C1-EXP. Go to step 4. 

4. Take the XOR of F1-SIGN and F2-SIGN to get 
C-SIGN. Go to step 5. 

5. Compute F1-HPF2-LO. Let the result be MI- 
HI. M1-LO. Goto step 6. 

6. Divide M1-HI.M1-LOby F2-HI. Let the quotient be 
M2-HI. Go to step 7. 

7. Do a true subtract of M2-HI from F1-LO. Let the 
result be M3-LO If C = 1 as a result of this sub- 
tract, then go to step 8, else decrement F1-HI and 
go to step 8. 

8. Divide F1-HI.M3-LO by F2-HI. Let the quotient be 
C1-HI and the remainder R1. Go to step 9. 

9. Divide R1 .0000 by F2-HI. Let the quotient be C1- 
LO. Go to step 1 0. 

10. If the MSB of C1-HI is 1 then go to step 11, else 
shift left C1-HI.C1-LO, decrease C1-EXP by 1 and 
go to step 1 1 . 

11. Round C1-HI.C1-LO to get C-HI.C-LO. go to step 
12. 

12. Pack C1-EXP.C-SIGN, C-HI and C-LO and return 
as the result. 

(v) ATOF. 

The processing steps in this case are as follows. 

1. Set M-SIGN, the mantissa sign to 0. 

Set M10-EXP, the implicit decimal exponent to 0. 
Set HI-INT to 0. 
Set LO-INT to 0. 
Go to step 2. 

2. Get a character from the input string. Let the char- 
acter be C. 

If C is a ' + ', then go to the start of step 2. 

If C is a ' - ', then set M-SIGN to FF and go to start 

of step 2. 

If C is a '.', then go to step 5. 

If C is none of the above, then go to step 3. 

3. Subtract 30 from C to get its integer value. Let this 
be I. Check and see if (HI-INT.LO-INT)*10 + 9 
can fit in 32 bits. If it can, then go to step 3.1 , else 
go to step 3.2. 

3.1. Multiply HI-INT.LO-INT by 10 and add I to the 
product. Store this sum back in HI-INT.LO-INT. Go 
to step 4. 

3.2. Increase M10-EXP by 1 and go to step 4. 

4. Get a character from the input string. Let the char- 
acter be C. 

If C is a '.', then go to step 5. 

If C is a 'E', then go to step 7. 

If C is the space character, then go to the start of 

step 4. 

If C is none of the above, then go to step 3. 



5. Get a character from the input string. Let the char- 
acter be C. 

If C is a 'E', then go to step 7. 

If C is the space character, then go to the start of 

step 5. 

If C is none of the above, then go to step 6. 

6. Subtract 30 from C to get its integer value. Let this 
be I. Check and see if (HI-INT.LO-INT)*10 + 9 
can fit in 32 bits. If it can, then go to step 6.1 , else 
go to step 5. 

6.1. Multiply HI-INT.LO-INT by 10 and add I to the 
product. Store this sum back in HI-INT.LO-INT. 
Decrement M10-EXP by 1. Go to step 5. 

7. Set SEXP, the exponent sign to be 0. Go to 
step 8. 

8. Get a character from the input string. Let the char- 
acter be C. 

If C is a ' + ', then go to start of step 8. 

If C is a ' - ', then set SEXP to be FF and go to the 

start of step 8. 

If C is none of the above, then go to step 9. 

9. Set M20-EXP, the explicit decimal exponent to 0. 
Go to step 10. 

1 0. Subtract 30 from C to get its integer value. Let this 
be I. Multiply M20-EXP by 10 and add I to the 
product. Store this sum back in M20-EXP. Go to 
step 1 1 . 

1 1 . Get a character from the input string. Let this be 
C. If C is the null character, then go to step 12, 
else go to step 1 0. 

12. Add M10-EXP and M20-EXP (with the proper sign 
as determined by SEXP) to get the 10's exponent 
M-EXP. Save in M-EXP the magnitude of the sum 
and in SEXP the sign of the sum. Go to step 13. 

13. Check and see if HI-INT.LO-INT is 0. If it is, then 
set the resulting floating point number, C, to zero 
and return. If it is not then go to step 14. 

14. Normalize HI-INT.LO-INT by left shifts such that 
the MSB is 1 . Let the number of left shifts needed 
to do this be L Set B1-EXP to 32-decimal - L Go 
to step 1 5. 

15. If SEXP is 0, then set P-HI.P-LO to the binary rep- 
resentation of 0.625, else set P-HI.P-LO to the bi- 
nary representation of 0.8. Go to step 1 6. 

16. Multiply HI-INT.LO-INT by P-HI.P-LO M-EXP 
times. After each multiplication, normalize the par- 
tial product if needed by left shifting. Accumulate 
the number of left shifts needed in B2-EXP. Let 
the final product be C-HI.C-LO. Go to step 17. 

17. Subtract B2-EXP from B1-EXP. Let the result be 
B-EXP. Go to step 18. 

1 8. If SEXP is 0, then multiply M-EXP by 4, else multi- 
ply M-EXP by -3. Let the result be B3-EXP. Go to 
step 1 9. 

19. Add B-EXP and B3-EXP. Let the result be C1- 
EXP. Add 126 to C1-EXP to restore the IEEE bias, 
getting C-EXP. Go to step 20. 

20. Round C-HI.C-LO. Go to step 21. 

21 . Pack C-EXP.M-SIGN, C-HI and C-LO and return, 
(vi) FTOA. 

The processing steps are as follows. 



1. Unpack the input FLP number. Let the unpacked 
number be represented by C-EXP.C-SIGN, C-HI 
and C-LO. Go to step 2. 

2. Subtract 126-decimal from C-EXP to remove the 
IEEE bias. Let the result be C1-EXP. Go to step 3. 

3. Multiply C1-EXP by the binary representation of 
log(2). Let the product be U-HI.U-LO Go to step 
4. 

4. Subtract 8 from U-HI.U-LO. Let the magitude of 
the integer part of the result be V and its sign 
VSIGN. Go to step 5. 

5. If VSIGN is 0, then set P-HI.P-LO to the binary 
representation of 0.8, else set P-HI.P-LO to the 
binary representation of 0.625. Go to step 6. 

6. Multiply C-HI.C-LO by P-HI.P-LO V times. Normal- 
ize the partial product after each multiplication, if 
needed, by left shifting. Accumulate any left shifts 
needed in B1-EXP. Let the final product be HI- 
INT.LO-INT. Go to step 7. 

7. Subtract B1-EXP from C1-EXP. Let the result be 
B2-EXP. Go to step 8. 

8. If VSIGN is 0, then multiply V by -3, else multiply 
it by 4. Let the result be B3-EXP. Go to step 9. 

9. Add B2-EXP and B3-EXP. Let the result be B4- 
EXP. Go to step 1 0. 

10. If B4-EXP is more than 32-decimal, then increase 

V and go to step 6, else go to step 1 1 . 

11. If B4-EXP is less than 28-decimal, then decrease 

V and go to step 6, else go to step 12. 

12. Subtract B4-EXP from 32. Let the result be B5- 
EXP. Goto step 13. 

13. Shift HI-INT.LO-INT right B5-EXP number of 
times. Go to step 14. 

14. Add 16-decimal to the address of the start of the 
decimal string. Output a null byte there. Go to step 
15. 

1 5. Divide V by 1 0-decimal. Let the quotient be Q and 
the remainder R. Add 30 to R and output it to the 
decimal string. Next add 30 to Q and output it to 
the decimal string. Go to step 16. 

16. If VSIGN is 0, then output ' + ' to the output string, 
else output ' - ' to the output string. Go to step 1 7. 

1 7. Output 'E' to the output string. Output '.' to the 
output string. Go to step 18. 

18. Divide C-HI.C-LO by 10-decimal 10 times. Let the 
remainder in each division be R. Add 30 to each R 
and output it to the output string. Go to step 19. 

19. If C-SIGN is 0, then output the space character to 
the output string, else output '-'to the output 
string. Then return to the calling program. 
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APPENDIX A 

A USER'S MANUAL FOR THE HPC 
FLOATING POINT PACKAGE 

The Single Precision Floating Point Package for the HPC 
implements the following functions. 

ARITHMETIC FUNCTIONS 

1 . FADD — Add two floating point numbers. 

2. FSUB — Subtract two floating point numbers. 

3. FMULT — Multiply two floating point numbers. 

4. FDIV — Divide two floating point numbers. 

FORMAT CONVERSION FUNCTIONS 

5. ATOF — Convert an ASCII string representing a decimal 
floating point number to a single precision floating point 
number. 

6. FTOA — Convert a single precision floating point number 
to an ASCII string that represents the decimal floating 
point value of the number. 

The entire package is in the form of a collection of subrou- 
tines and is contained in the following files. 

1. FERR.MAC 

2. FNACHK.MAC 

3. FZCHK.MAC 

4. FUNPAK.MAC 

5. FPAK.MAC 

6. FPTRAP.MAC 

7. ROUND.MAC 

8. BFMULMAC 

9. ISIOK.MAC 

10. MUL1 O.MAC 

1 1 . ATOF.MAC 

12. FTOA.MAC 

13. FADD.MAC 

14. FMULT.MAC 

15. FDIV.MAC 

The first 7 files are general utility routines that are used by 
all the Arithmetic and Format Conversion subroutines. The 
next 3 files, BFMULMAC, ISIOK.MAC and MUL1 O.MAC are 
used only by the Format Conversion subroutines, ATOF and 
FTOA. Depending on the functions being used in the user 
program, only the necessary files need be included. 

INTERFACE WITH USER PROGRAMS 

1. All the Arithmetic routines expect the input to be in the 
IEEE Single Precision format. This format requires 2 
words for the storage of each floating point number. If the 
required arithmetic operation is FlopF2, where op is +, 
— , * or /, then the routines expect that F1 is available in 
registers K and A on entry, with the high half in K. Also, 
the two words of F2 are expected to be on the stack. If 
SP is the stack pointer on entry into one of the Arithmetic 
function subroutines, then the high word of F2 should be 
at W(SP-4) and the low word at W(SP-6). The result of the 
Arithmetic operation is returned in IEEE format in regis- 
ters K and A, with the high word in K. 



2. The Format Conversion subroutine ATOF expects that on 
entry, B contains the address of the ASCII string repre- 
senting the decimal floating point number. This string 
must be of the form 

Siiiii.ffffffEsNND 
where 

S is an optional sign for the mantissa. Thus S can be ' + ', 
'-'or not present at all. 

iiiii is the optional integer part of the mantissa. If it is pres- 
ent, it can be of any length, must contain only the charac- 
ters '0' through '9' and must not contain any embedded 
blanks. 

. is the optional decimal point. It need not be present if 
the number has no fractional part, 
ffffff is the optional fractional part of the mantissa, ffffff, if 
it is present must consist of a sequence of digits '0' 
through '9'. It can be of any length. Note that either iiiii, 
the integer part or .ffffff the fractional part must be pres- 
ent. 

E is the required exponent start symbol, 
s is the optional sign of the exponent. If it is present, it 
must be ' + ' or ' — '. 

NN is the exponent and consists of at most two decimal 
digits. It is required to be present. 
D is the null byte <00> and must be present to terminate 
the string. 

The floating point number represented by the above 
string is returned by ATOF in IEEE format in registers K 
and A. 

3. The format conversion routine FTOA expects the floating 
point number input to be in registers K and A in the IEEE 
format. Register B is expected to contain the starting ad- 
dress of a 17 byte portion of memory where the output 
string will be stored. 

4. Three global symbols need to be defined in the user pro- 
gram before assembling the user program and any includ- 
ed floating point package files. These symbols are: 

(i) TMP1 which must be set to 0. The package uses 

W(TMP1) for temporary storage, 
(ii) FPERWD which must be set to an address in the base 

page. The package signals floating point exceptions 

using W(FPERWD). This is described below, 
(iii) FPTRAP which must be set to the address of the start 

of a user floating point exception handler. Again this is 

described below. 

FLOATING POINT EXCEPTS 

The package maintains a history of floating point exceptions 
in the 4 least significant bits of the word W(FPERWD). The 
value of the symbol FPERWD should be defined by the user 
program, and should be an address in the base page. This 
word should also be cleared by the user program before 
calling any floating point routine. The word is never cleared 
by the floating point package, and the user program can 
examine this word to determine the type of exceptions that 
may have occurred during the course of a computation. 



The following 4 types of error can occur in the course of a 
floating point computation. 

1. Invalid Operand. This happens if one of the input num- 
bers for an Arithmetic routine or the input for FTOA is not 
a valid floating point number. An invalid floating point 
number (or NaN) can be created either by an overflow in 
a previous computation step, or if the ASCII decimal float- 
ing point number input to ATOF is too large to be repre- 
sented in the IEEE format. The result, if one of the inputs 
is a NaN is always set to a NaN. 

2. Overflow. This happens if the result of a computation is 
too large to be represented within the exponent range 
available. Overflow can occur in any of the arithmetic rou- 
tines or ATOF. On overflow, the result is set to a repre- 
sentation called NaN. An NaN is considered an illegal 
operand in all successive steps. 

3. Underflow. This occurs if the result of a computation is 
too small to be represented with the precision and expo- 



nent range available. On underflow, the result is set to 
zero. 
4. Divide-by-zero. This error occurs if F2 is zero when com- 
puting F1/F2. The result is set to an NaN. 
Each of the above errors results in a bit being set in 
W(FPERWD). This is done as follows: 
Bit 0— Set on Overflow. 
Bit 1— Set on Underflow. 
Bit 2— Set on Illegal Operand. 
Bit 3 — Set on Divide-by-zero. 

One further action is taken when a floating point exception 
occurs. After the result has been set to the appropriate val- 
ue, and the corresponding bit in W(FPERWD) set, the pack- 
age does a subroutine call to address FPTRAP. The user 
can provide any exception handler at this address. The file 
FPTRAP.MAC contains the simplest possible user excep- 
tion handler. It does nothing, but merely returns back to the 
calling program. 
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FLP 



1 






.TITLE FLP 


2 




LISTER: 




3 


0071 




.LIST 071 


4 


FOOO 




. = OFOOO 


5 


0002 




FPERWD = W(2) 


6 


0000 




TMP1 = W(0) 
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FLP 

THE FLP ROUTINES 



.FORM 'THE FLP ROUTINES' 



The code listed in this App Note is available on Dial-A-Helper. 

Dial-A-Helper is a service provided by the Microcontroller Applications Group. The Dial-A-Helper system provides access to 
an automated information storage and retrieval system that may be accessed over standard dial-up telephone lines 24 hours 
a day. The system capabilities include a MESSAGE SECTION (electronic mail) for communicating to and from the Microcon- 
troller Applications Group and a FILE SECTION mode that can be used to search out and retrieve application data about 
NSC Microcontrollers. The minimum system requirement is a dumb terminal, 300 or 1200 baud modem, and a telephone. 
With a communications package and a PC, the code detailed in this App Note can be down loaded from the FILE SECTION 
to disk for later use. The Dial-A-Helper telephone lines are: 

Modem (408) 739-1162 

Voice (408)721-5582 

For Additional Information, Please Contact Factory 



NATIONAL SEMICONDUCTOR 


CORPORATION PAGE 


: 3 


HPC CROSS 


ASSEMBLER, REV 


:C,30 JUL 


86 




FLP 












FERR 


MAC 










8 








.FORM ' 


FERR. MAC 


9 








.INCLD 


FERR. MAC 


1 






; EXCEPTION HANDLING. 




2 






; DIVIDE BY ZERO. 




3 






DIVBYO : 






4 


FOOO 


820802FA 




OR FPERWD, 08 


; SET THE DIVIDE BY BIT. 


5 


F004 


00 




CLR A 




6 


F005 


B17F80 




LD K, 07F80 




7 


F008 


3093 




JSR FPTRAP 




8 


FOOR 


3FCC 




POP B 




9 


FOOC 


3FCE 




POP X 




10 


FOOE 


3C 




RET 




11 






; ILLEGAL OPERAND - ONE 


OF Fl OR F2 IS A NAN. 


12 






FNAN: 






13 


FOOF 


820402FA 




OR FPERWD, 04 


; SET THE ILLEGAL OPERAND BIT. 


14 


F013 


00 




CLR A 




15 


F014 


B17F80 




LD K, 07F80 


; RETURN NAN IN K AND A. 


16 


F017 


3084 




JSR FPTRAP 


; GO TO USER TRAP ROUTINE. 


17 


F019 


3FCC 




POP B 




18 


F01B 


3FCE 




POP X 




19 


F01D 


3C 




RET 




20 






; EXPONENT UNDERFLOW. 




21 






UNDFL : 






22 


F01E 


820202FA 




OR FPERWD, 02 


; SET THE EXPONENT UNDERFLOW BIT. 


23 


F022 


00 




CLR A 




24 


F023 


ACC8CA 




LD K, A 




25 


F026 


3075 




JSR FPTRAP 




26 


F028 


3FC4 




POP SP 




27 


F02A 


3FCC 




POP B 




28 


F02C 


3FCE 




POP X 




29 


F02E 


3C 




RET 




30 






; EXPONENT OVERFLOW. 




31 






OVRFL : 






32 


F02F 


820102FA 




OR FPERWD, 01 


; SET THE EXPONENT OVERFLOW BIT. 


33 


F033 


00 




CLR A 




34 


F034 


B17F80 




LD K, 07F80 




35 


F037 


3064 




JSR FPTRAP 




36 


F039 


3FC4 




POP SP 




37 


F03B 


3FCC 




POP B 




38 


F03D 


3FCE 




POP X 




39 


F03F 


3C 




RET 




40 






; 






41 








.END 
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FLP 
FNACHK.MAC 



PAGE: 



10 
11 

1 
2 
3 
4 
5 



9 
10 
11 
12 

13 F040 AECA 

14 F042 E7 

15 F043 BDFEFF 

16 F046 45 

17 F047 D7 

18 F048 03 

19 F049 AECA 

20 F04B 3C 
21 

22 F04C D7 

23 F04D 02 

24 F04E AECA 

25 F050 3C 
26 

27 



.FORM 'FNACHK.MAC 
.INCLD FNACHK.MAC 
.TITLE FNACHK 
.LOCAL 

SUBROUTINE TO CHECK IF A SP FLOATING POINT NUMBER STORED IN THE 
IEEE FLOATING POINT FORMAT IN REGS. K AND A IS NAN. 



RETURNS IN 
RETURNS 1 IN 



IF THE NUMBER IS NOT A NAN. 
IF THE NUMBER IS A NAN. 



PRESERVES REGS. K, A, X AND B. DESTROYS C. 



FNACHK : 


X A, K 




SHL A 




IFGT A.OFEFF 




JP $ISNAN 




RRC A 




RESET C 




X A, K 




RET 


$ISNAN: 






RRC A 




SET C 




X A, K 




RET 



.END 
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FZCHK.MAC 



PAGE: 



12 






13 






1 






2 






3 






4 






5 






6 






7 






8 






9 






10 






11 






12 


F051 


AECA 


13 


F053 


E7 


14 


F054 


9DFF 


15 


F056 


45 


16 


F057 


D7 


17 


F058 


02 


18 


F059 


AECA 


19 


F05B 


3C 


20 






21 


F05C 


D7 


22 


F05D 


03 


23 


F05E 


AECA 


24 


F060 


3C 


25 






26 







.FORM 'FZCHK.MAC 
.INCLD FZCHK.MAC 
.TITLE FZCHK 
.LOCAL 

SUBROUTINE THAT CHECKS IF A SP FLOATING POINT NUMBER STORED 
IN THE IEEE FORMAT IN REGS K AND A IS ZERO. 

RETURNS IN C IF THE NUMBER IS NOT ZERO. 
RETURNS 1 IN C IF THE NUMBER IS ZERO. 
SAVES REGS. K, A, X, AND B BUT DESTROYS C. 



FZCHK: 



$ANOTO : 



X A, K 


SHL A 


IFGT A, OFF 


JP }ANOTO 


RRC A 


SET C 


X A, K 


RET 


RRC A 


RESET C 


X A, K 


RET 



10 
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FZCHK 








FUNPAK.MAC 






14 






.FORM 


FUNPAK.MAC 


15 






.INCLD 


FUNPAK.MAC 


1 






.TITLE 


FUNPAK 


2 
3 
4 






.LOCAL 








SUBROUTINE TO UNPACK 


A SP FLOATING POINT NUMBER STORED IN THE 


5 






IEEE FORMAT IN REGS. 


K AND A. THE UNPACKED FORMAT OCCUPIES 3 


6 
7 
8 






WORDS AND IS ORGANIZED AS FOLLOWS: 






increasing addrs | 


<- X on exit 


9 








|eeeeeeeessssssss| fexp-fsign 


10 








JmmmmmmmmmmmmmmmmJ fhi 


11 








MMMMMMMMOOOOOOOO FLO <- X on entry 


12 










13 










14 






EEEEEEEE - 8 BIT EXPONENT IN EXCESS-127 FORMAT 


15 






SSSSSSSS - SIGN BIT 


< 00 -> +, FF -> -> 


16 






M ... M - 24 BITS OF 


MANTISSA. NOTE THAT IMPLIED 1 IS PRESENT HERE. 


17 










18 






ON ENTRY TO THE SUBROUTINE X SHOULD POINT TO FLO. ON EXIT, X POINTS 


19 






TO THE WORD AFTER FSIGN. 


20 






REGS. K, A AND B ARE 


DESTROYED BY THIS SUBROUTINE. 


21 










22 




J 


UNPAK : 




23 


F061 


ABCC 


ST A,B 


; SAVE A IN B. 


24 


F063 


00 


CLR A 




25 


F064 


Dl 


X A, M(X+) 


; ZERO LOW BYTE OF FLO. 


26 


F065 


88CC 


LD A, L(B) 




27 


F067 


Dl 


X A, M(X+) 


; MOVE LOW BYTE OF F-RO INTO HIGH BYTE OF FLO. 


28 


F068 


88CD 


LD A, H(B) 




29 


F06A 


Dl 


X A, M(X+) 


; MOVE MID BYTE OF MANT INTO LOW BYTE OF FHI. 


30 


F06B 


A8CA 


LD A, K 




31 


F06D 


96C80F 


SET A. 7 


; SET IMPLIED 1 IN MANTISSA 


32 


F070 


Dl 


X A, M(X+) 


; MOVE HIGH BYTE OF MANT INTO HIGH BYTE OF FHI. 


33 


F071 


A8CA 


LD A, K 




34 


F073 


E7 


SHL A 


; SIGN BIT TO CARRY. 


35 


F074 


B9FF00 


AND A, OFFOO 


; ZERO SIGN. 


36 


F077 


07 


IF C 




37 


F078 


9AFF 


OR A, OFF 


; PUT SIGN BACK IF -. 


38 


F07A 


Fl 


X A, W(X+) 


; SAVE FEXP-FSIGN. 


39 


F07B 


3C 


RET 




40 










41 






.END 





11 
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FUNPAK 








FPAK 


MAC 








16 






.FORM 


FPAK. MAC 


17 






.INCLD 


FPAK. MAC 


1 






.TITLE 


FPAK 


2 
3 

4 






.LOCAL 










SUBROUTINE TO PACK A 


SP FLOATING POINT NUMBER STORED IN THE 


5 








3 WORD FEXP-FSIGN/FHI/FLO FORMAT INTO THE IEEE FORMAT IN REGS. 


6 
7 
8 








K AND A. 










ON ENTRY TO THE SUBROUTINE, X POINTS TO FLO. ON EXIT, X POINTS 


9 








TO THE WORD AFTER FSIGN. 


10 












11 








REGS. K, A AND B ARE 


DESTROYED. 


12 












13 






FPAK: 




14 


F07C Dl 




X A, M(X+) 




GET RID OF ZERO LOW BYTE OF FLO. 


15 


F07D Dl 




X A, M(X+) 




GET HIGH BYTE OF FLO. 


16 


F07E ABCA 




ST A, K 




STORE IT IN K. 


17 


F080 Dl 




X A, M(X+) 




GET LOW BYTE OF FHI . 


18 


F081 3B 




SWAP A 




19 


F082 3B 




SWAP A 




20 


F083 B9FF00 




AND A, OFFOO 


; SHIFT LEFT 8 TIMES. 


21 


F086 A0C8CAFA 




OR K, A 


; LOW WORD OF RESULT IS NOW IN K. 


22 






; 




23 


F08A Dl 




X A, M(X+) 




GET HIGH BYTE OF FHI. 


24 


F08B 96C81F 




RESET A. 7 




ZERO IMPLIED MSB 1 IN MANT. 


25 


F08E ABCC 




ST A,B 




SAVE IN REG. B. 


26 


F090 D4 




LD A, M(X) 




GET SIGN BYTE FROM ASIGN. 


27 


F091 C7 




SHR A 




MOVE 1 SIGN BIT INTO CARRY. 


28 


F092 FO 




LD A, W(X+) 




GET FEXP-FSIGN. 


29 


F093 B9FF00 




AND A, OFFOO 




ZERO SIGN. 


30 


F096 D7 




RRC A 




MOVE RIGHT 1 BIT. SIGN BIT FROM C 


31 










ENTERS INTO THE MSB. 


32 


F097 96CCFA 




OR A, B 




GET MANT BITS IN FROM B. 


33 


F09A AECA 




X A, K 




SWAP A AND K 


34 


F09C 3C 




RET 




35 






; 




36 






.END 
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HPC CROSS ASSEMBLER 


REV 
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FPAK 










FPTRAP.MAC 








18 






.FORM 


FPTRAP.MAC 


19 






.INCLD 


FPTRAP.MAC 


1 






.TITLE 


FPTRAP 


2 






; USER SUPPLIED FP TRAP ROUTINE. 


3 






FPTRAP : 




4 
5 
6 


F09D 3C 




RET 








.END 





12 
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FPTRAP 










ROUND. MAC 










20 








.FORM 'ROUND. MAC 


21 








.INCLD ROUND. MAC 


1 








.TITLE SROUND 


2 
3 
4 








.LOCAL 








THIS 


SUBROUTINE IS USEI 


TO ROUND THE 32 BIT MANTISSA OBTAINED 


5 
6 
7 






IN THE FLOATING POINT CALCULATIONS TO 24 BITS. 






THE UNPACKED FLOATING POINT NUMBER SHOULD BE STORED IN 


8 






CONSECUTIVE WORDS OF MEMORY. ON ENTRY, X SHOULD CONTAIN 


9 






THE ADDRESS OF C-HI. C- 


EXP. C-SIGN IS AT W(X+2) AND 


10 






C-LO 


IS AT W(X-2) . 




11 












12 






ON EXIT X HAS THE ADDRESS OF C-EXP. C-SIGN. 


13 




c 


ROUND : 






14 


F09E 


F2 




LD A, W(X-) 


REMEMBER X POINTS TO C-HI. 


15 


F09F 


F4 




LDA, W(X) 


LOAD C-LO. 


16 


FOAO 


96C817 




IF A. 7 


IF BIT 25 OF MANTISSA IS 1, 


17 


F0A3 


43 




JP $RNDUP 


THEN NEED TO INCREASE MANTISSA. 


18 


F0A4 


FO 




LD A, W(X+) 




19 


F0A5 


FO 




LD A, W(X+) 


X NOW POINTS TO C-EXP. C-SIGN. 


20 


F0A6 


5F 




JP $EXIT 


DONE, SO GET OUT. 


21 






INCREASE MANTISSA. 




22 




1 


RNDUP : 






23 


F0A7 


B80100 




ADD A, 0100 




24 


FOAA 


Fl 




X A, W(X+) 


INCREASE LOW BYTE BY 1. 


25 


FOAB 


07 




IF C 


IF THERE IS A CARRY, 


26 


FOAC 


42 




JP JHIUP 


THEN NEED TO INCREASE C-HI. 


27 


FOAD 


FO 




LD A, W(X+) 


X NOW POINTS TO C-EXP. C-SIGN. 


28 


FOAE 


57 




JP $EXIT 


DONE, SO GET OUT. 


29 






MANTISSA INCREASE PROPAGATING TO HIGH WORD. 


30 




1 


HIUP: 






31 


FOAF 


F4 




LD A, W(X) 




32 


FOBO 
F0B1 
F0B2 


B8 
00 
01 




.BYTE 0B8,00,01 


DO ADD A, 01 BUT WITH WORD CARRY!! 


33 


F0B3 


07 




IF C 


IF THERE IS A CARRY, 


34 


F0B4 


42 




JP $EXIN2 


THEN NEED TO INCREASE EXPONENT. 


35 


F0B5 


Fl 




X A,W(X+) 




36 


F0B6 


4F 




JP $EXIT 


GET OUT. 


37 






ROUND 


UP LEADS TO EXPOI* 


ENT INCREASE. 


38 




1 


EXIN2 : 






39 


F0B7 


D7 




RRC A 


CARRY->MSB, LSB->CARRY. 


40 


F0B8 


F3 




X A,W(X-) 




41 


F0B9 


F4 




LD A,W(X) 


LOW WORD IS NOW IN A. 


42 


FOBA 


D7 




RRC A 




43 


FOBB 


Fl 




X A, W(X+) 




44 


FOBC 


FO 




LD A, W(X+) 


X NOW POINTS TO C-EXP. CSIGN. 


45 


FOBD 


F4 




LD A, W(X) 




46 


FOBE 


B80100 




ADD A, 0100 




47 


FOCI 


07 




IF C 
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SROUND 



OR A, OFFOO ; MAKE IT A NAN. 
ST A, W(X) 



ROUND. MAC 




48 F0C2 


BAFFOO 


49 F0C5 


F6 


50 




51 




52 F0C6 


3C 


53 




54 





$EXII: 

RET 



.END 



14 
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BFMUL.MAC 



PAGE: 



22 
23 
1 
2 
3 
4 
5 



9 






10 






11 






12 






13 






14 






15 






16 






17 






18 






19 






20 






21 


F0C7 


AFCE 


22 


F0C9 


AFC8 


23 


FOCB 


AFC A 


24 


FOCD 


ABC A 


25 


FOCF 


A6FFF6C4FE 


26 


F0D4 


3FCA 


27 


F0D6 


AFCE 


28 


F0D8 


AFC8 


29 


FODA 


A8CA 


30 


FODC 


A6FFF2C4FE 


31 


F0E1 


3FC8 


32 






33 


F0E3 


3FCA 


34 


F0E5 


96CEF8 


35 


F0E8 


07 


36 


F0E9 


A9CA 


37 


FOEB 


3FCE 


38 


FOED 


AFC A 


39 


FOEF 


AFC8 


40 


F0F1 


A8CE 


41 


F0F3 


A6FFF6C4FE 


42 


F0F8 


3FC8 


43 


FOFA 


3FCA 


44 


FOFC 


96CEF8 


45 


FOFF 


07 


46 


F100 


A9CA 


47 


F102 


3FCE 


48 


F104 


3C 


49 







.FORM 'BFMUL.MAC 
.INCLD BFMUL.MAC 
.TITLE BFMUL 

THIS SUBROUTINE IS USED TO MULTIPLY TWO 32 BIT FIXED POINT FRACTIONS. 
THE ASSUMED BINARY POINT IS TO THE IMMEDIATE LEFT OF THE MSB. 

THE FIRST FRACTION IS STORED IN REGS K AND A, WITH THE MORE 
SIGNIFICANT WORD BEING IN K. 

THE SECOND FRACTION IS STORED ON THE STACK. THE MORE SIGNIFICANT 
WORD IS AT W(SP-4) AND THE LOWER SIGNIFICANT WORD 
IS IN THE WORD BELOW II. 

THE 32 BIT PRODUCT IS LEFT IN REGS. K AND A, WITH THE MORE 
SIGNIFICANT WORD BEING IN K. 

IMPORTANT NOTE : THE FRACTIONS ARE ASSUMED TO BE UNSIGNED. 

REGS. B AND X ARE UNCHANGED. 



BFMUL : 



PUSH X 




PUSH A 




PUSH K 




LD A, K 




MULT A, 


W(SP-OA) 


POP K 




PUSH X 




PUSH A 




LD A, K 




MULT A, 


W(SP-OE) 


POP A 




POP K 




ADD A, X 


IF C 




INC K 




POP X 




PUSH K 




PUSH A 




LD A, X 




MULT A, 


W(SP-OA) 


POP A 




POP K 




ADD A, X 


IF C 




INC K 




POP X 




RET 





SAVE X. 

SAVE Fl-LO 

SAVE Fl-HI. 

MOVE Fl-HI TO A. 

MULTIPLY Fl-HI BY F2-HI . 

GET Fl-HI. 

SAVE PR-HI. 

SAVE PR-LO. 

MOVE Fl-HI TO A. 

MULTIPLY Fl-HI BY F2-L0. 

GET PR-LO SAVED. NOTE THAT THE 

LO WORD OF THIS PRODUCT IS DISCARDED. 

GET PR-HI SAVED. 

ADD TO PR-LO THE HI WORD OF THIS PRODUCT. 

ON CARRY, 

PROPAGATE THRU TO PR-HI. 

GET Fl-LO. 

SAVE PR-HI. 

SAVE PR-LO. 

MOVE Fl-LO TO A. 

MULTIPLY BY F2-HI . 

GET PR-LO SAVED. 

GET PR-HI SAVED. 

ADD TO PR-LO THE HI-WORD OF THIS PRODUCT. 

PROPAGATE ANY CARRY TO PR-HI. 
RESTORE X. 
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BFMUL 
BFMUL. MAC 



PAGE: 



12 



50 



.END 



NATIONAL SEMICONDUCTOR CORPORATION 
HPC CROSS ASSEMBLER, REV :C, 30 JUL 86 
BFMUL 
ISIOK.MAC 



PAGE: 



13 



24 




25 




1 




2 




3 




4 




5 




6 




7 




8 




9 




10 


I 


11 F105 


02 


12 F106 


861999CAFC 


13 F10B 


47 


14 F10C 


861999CAFD 


15 Fill 


03 


16 F112 


3C 


17 F113 


BD9998 f) 


18 F116 


03 


19 F117 


3C 


20 




21 





.FORM 'ISIOK.MAC 
.INCLD ISIOK.MAC 
.TITLE ISIOK 
.LOCAL 

THIS SUBROUTINE IS USED TO DETERMINE IF ANOTHER DECIMAL DIGIT CAN 
BE ACCUMULATED IN THE 32 BIT INTEGER STORED IN REGS. K AND A. 
THE MORE SIGNIFICANT WORD IS IN K. 

SETS THE CARRY TO 1 IF IT CAN BE ACCUMULATED ; RESETS THE CARRY 
OTHERWISE. PRESERVES ALL REGS. 



ISIOK: 



$CHKOT 



SET C 

IFEQ K, 01999 

JP ftCHKOT 

IFGT K, 01999 

RESET C 

RET 

IFGT A, 09998 

RESET C 

RET 

.END 
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26 
27 
1 
2 
3 
4 
5 



.FORM 'MUL10.MAC 
.INCLD MUL10.MAC 
.TITLE MUL10 
.LOCAL 

THIS SUBROUTINE MULTIPLIES THE 32 BIT INTEGER STORED IN REGS K AND A 
BY 10-DECIMAL AND ADDS TO IT THE INTEGER STORED IN X. 
THE RESULT IS RETURNED IN K AND A. 
REGS. B AND X ARE NOT CHANGED. 



MUL10 : 



10 


F118 


AFCE 


11 


F11A 


AFC8 


12 


F11C 


A8CA 


13 


F11E 


9E0A 


14 


F120 


AFC8 


15 


F122 


A6FFFCC4A8 


16 


F127 


9E0A 


17 


F129 


3FCA 


18 


F12B 


A0CECAF8 


19 


F12F 


3FCE 


20 


F131 


3FCE 


21 


F133 


96CEF8 


22 


F136 


07 


23 


F137 


A9CA 


24 


F139 


3C 


25 






26 






28 







PUSH X 
PUSH A 

LD A 
MULT 
PUSH 
LD A 



K 



OA 



W(SP-4) 
MULT A, OA 
POP K 



ADD K, 
POP X 
POP X 
ADD A, 
IF C 
INC K 
RET 



X 



SAVE INTEGER. 
SAVE LONG INT-LO. 

MULT LONG INT-HI BY 10. 
SAVE LOW WORD OF PRODUCT. 
GET LONG INT-LO. 

GET LO WORD OF LAST PRODUCT. 

ADD TO IT HI WORD OF THIS PRODUCT. 

GET RID OF GARBAGE. 

GET INTEGER TO BE ADDED. 



.END 
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.FORM 'ATOF. MAC 


30 






.INCLD ATOF. MAC 


1 






.TITLE ATOF 


2 
3 

4 






.LOCAL 






THIS SUBROUTINE CONVERTS A DECIMAL FLOATING POINT STRING TO 


5 






AN IEEE FORMAT SINGLE PRECISION FLOATING POINT NUMBER. THE 


6 






INPUT DECIMAL STRING IS ASSUMED TO BE OF THE FORM 


7 






SMMMMMMM . FFFFFEDNN 


8 






WHERE S IS THE SIGN OF THE DECIMAL MANTISSA, 


9 






M...M IS THE INTEGER PART OF THE MANTISSA, 


10 






F...F IS THE FRACTIONAL PART OF THE MANTISSA, 


11 






D IS THE SIGN OF THE DECIMAL EXPONENT, 


12 






AND NNN IS THE DECIMAL EXPONENT. 


13 








14 






ON ENTRY, B SHOULD POINT TO THE ADDRESS OF THE ASCII 


15 






STRING HOLDING THE DECIMAL FLOATING POINT NUMBER. THIS STRING 


16 






MUST BE TERMINATED BY A NULL BYTE. 


17 








18 






THE BINARY FLOATING POINT NUMBER IS RETURNED IN 


19 






REGS. K AND A. 


20 








21 






REGS. B AND X ARE LEFT UNCHANGED. 


22 








23 








24 








25 








26 








27 




ATOF: 


28 


F13A 


AFCE 


PUSH X 


29 


F13C 


AFCC 


PUSH B 


30 


F13E 


00 


CLR A 


ZERO A. 


31 


F13F 


AFC8 


PUSH A 


STORAGE FOR MANTISSA SIGN. 


32 


F141 


AFC8 


PUSH A 


STORAGE FOR IMPLICIT 10'S EXPONENT. 


33 


F143 


AFC8 


PUSH A 


STORAGE FOR HI-INT. 


34 


F145 


AFC8 


PUSH A 


STORAGE FOR LO-INI. 


35 








36 






DECIMAL STRING MUST START WITH A ' + ', '-', '.' OR A DIGIT. 


37 






RESULTS ARE UNPREDICTABLE IF IT DOES NOT. 


38 






THE •+' MEANS THAT THE MANTISSA IS POSITIVE. IT CAN BE OMITTED. 


39 






THE '-' MEANS THAT THE MANTISSA IS NEGATIVE. 


40 






THE '.' MEANS THAT THE MANTISSA HAS NO INTEGER PART. 


41 








42 




J 


L00P1: 


43 


F147 


CO 


LDS A, M(B+) 


44 


F148 


40 


NOP 


GET THE CHARACTER. 


45 


F149 


9C2B 


IFEQ A, ■+' 


IF IT IS A '+' , 


46 


F14B 


64 


JP $L00P1 


DO NOTHING, BUT GET 1 MORE. 


47 


F14C 


9C2D 


IFEQ A, ■-' 


IF IT IS A ■-' , 


48 


F14E 


45 


JP $MSIGN 


GO AND CHANGE THE MANTISSA SIGN. 


49 


F14F 


9C2E 


IFEQ A, ' .' 


IF IT IS A ■ . ■ , 
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50 F151 


9438 


51 




52 F153 


48 


53 




54 F154 


90FF 


55 F156 


A6FFF8C4AB 


56 F15B 


74 


57 




58 




59 




60 




61 F15C 


02 


62 F15D 


8230C8EB 


63 F161 


ACC8CE 


64 F164 


3FCA 


65 F166 


3FC8 


66 F168 


3463 


67 F16A 


07 


68 F16B 


4B 


69 




70 




71 F16C 


3FCE 


72 F16E 


A9CE 


73 F170 


AFCE 


74 F172 


AFC8 


75 F174 


AFC A 


76 F176 


46 


77 




78 




79 




80 F177 


345F 


81 F179 


AFC8 


82 F17B 


AFC A 


83 




84 




85 F17D 


CO 


86 F17E 


40 


87 F17F 


9C2E 


88 F181 


49 


89 F182 


9C45 


90 F184 


9434 


91 F186 


9C20 


92 F188 


6B 


93 




94 F189 


952D 


95 




96 




97 




98 




99 F18B 


CO 


100 F18C 


40 



$MSIGN: 



JMP $FRCOL 

JP $INCOL 

LD A, OFF 

ST A, W(SP-08) 

JP $L00P1 



GO AND COLLECT THE FRACTION PART. 

GET HERE MEANS IT IS A DIGIT. 

SO GO AND COLLECT THE INTEGER PARI. 



CHANGE MANTISSA SIGN TO NEG. 
GO BACK FOR SOME MORE. 



$INCOL: 
; GET HERE MEANS COLLECTING INTEGER PART OF MANTISSA. 



SET C 

SUBC A, '0 
LD X, A 
POP K 
POP A 
JSR ISIOK 
IF C 
JP $ACCM 



POP X 
INC X 
PUSH X 
PUSH A 
PUSH K 
JP $ISNXT 



CONVERT DIGIT FROM ASCII TO INTEGER. 

MOVE INTEGER TO X. 

GET HI-INT COLLECTED SO FAR. 

GET LO-INT COLLECTED SO FAR. 

CHECK IF THE DIGIT CAN BE ACCUMULATED. 

LOOK AT C. 

YES, IT CAN BE SO GO DO IT. 

GET HERE MEANS CAN ACCUMULATE ANY MORE. 

SO INCREASE THE IMPLICIT 10 ' S EXPONENT. 

GET IMPLICIT 10 'S EXPONENT COLLECTED 

SO FAR AND INCREMENT IT. 

SAVE IT BACK. 

SAVE LO-INT. 

SAVE HI-INT. 



$ACCM: 
; GET HERE MEANS THE PRESENT DIGIT CAN BE ACCUMULATED. 



JSR MUL10 

PUSH A 

PUSH K 
$ISNXT: 
; PROCESS THE NEXT CHARACTER 

LDS A, M(B+) 

NOP 

IFEQ A, ' . ■ 

JP $FRCOL 

IFEQ A, 'E' 

JMP JEXCOL 

IFEQ A, ' ■ 

JP $ISNXT 



MULTIPLY BY 10 AND ADD DIGIT. 
SAVE LO-INT. 
SAVE HI-INT. 



IF IT IS A ' . ' 

GO COLLECT FRACTION PART. 

IF IT IS 'E' , 

GO COLLECT EXPONENT PART. 

IF IT IS A SPACE, 

GO GET SOME MORE. 

GET HERE MEANS IT IS A DIGIT. 



JMP SINCOL 



$FRCOL: 
; GET HERE MEANS COLLECT THE FRACTIONAL PART OF THE MANTISSA. 

LDS A, M(B+) 
NOP 
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101 F18D 


9C45 


IFEQ A, 'E' 




IF IT IS A 'E' , 


102 F18F 


9429 


JMP $EXCOL 




GO COLLECT EXPONENT. 


103 F191 


9C20 


IFEQ A, ' ' 




IF IT IS SPACE, 


104 F193 


68 


JP $FRCOL 




GO GET SOME MORE. 


105 








GET HERE MEANS IT IS A DIGIT. 


106 F194 


D2 


SET C 




107 F195 


8230C8EB 


SUBC A, '0' 




GET INTEGER FROM DIGIT. 


108 F199 


ACC8CE 


LD X, A 




SAVE IT IN A. 


109 F19C 


3FCA 


POP K 




GET HI-INT. 


110 F19E 


EFC8 


POP A 




GET LO-INT. 


111 F1A0 


349B 


JSR ISIOK 




CHECK IF IT CAN BE ACCUMULATED. 


112 F1A2 


07 


IF C 




113 F1A3 


45 


JP $ACCF 


; YES, SO GO DO II. 


114 






; GET HERE MEANS CAN'T COLLECT MORE DIGITS. 


115 F1A4 


AFC8 


PUSH A 




116 F1A6 


AFC A 


PUSH K 




117 F1A8 


7D 


JP $FRCOL 


; SO JUST IGNORE IT. 


118 




; 




119 




$ACCF: 




120 




; ACCUMULATE THE FRACTIONAL DIGIT. 


121 F1A9 


3491 


JSR MUL10 




MULTIPLY BY 10 AND ADD DIGIT. 


122 F1AB 


3FCE 


POP X 




GET IMPLICIT 10' S EXPONENT COLLECTED SO FAR, 


123 F1AD 


86FFFFCEF8 


ADD X, OFFFF 




AND DECREMENT IT BY 1. 


124 F1B2 


AFCE 


PUSH X 




SAVE IT BACK. 


125 F1B4 


AFC8 


PUSH A 




SAVE LO-INT. 


126 F1B6 


AFC A 


PUSH K 




SAVE HI-INT. 


127 F1B8 


952D 


JMP $FRCOL 




GO GET SOME MORE. 


128 




; 




128 




JEXCOL : 




130 




; GET HERE MEANS THE 


EXPLICIT 10 'S EXPONENT NEEDS TO BE 


131 




; COLLECTED FROM THE 


STRING. 


132 F1BA 


03 


RESET C 


; MAKE EXPONENT SIGN POST. 


133 




$EXCHR : 




134 F1BB 


CO 


LDS A, M(B+) 




135 F1BC 


40 


MOP 




136 F1BD 


9C2B 


IFEQ A, '+' 




IF IT IS A '+' , 


137 F1BF 


64 


JP $EXCHR 




GET SOME MORE. 


138 F1C0 


9C2D 


IFEQ A, ■-' , 




IF IT IS A ■-' , 


139 F1C2 


44 


JP $ESIGN 




GO FIX EXPONENT SIGN. 


140 F1C3 


9C20 


IFEQ A, ■ ' 




IF IT IS SPACE, 


141 F1C5 


6A 


JP $EXCHR 




GO GET SOME MORE. 


142 








GET HERE MEANS IT IS A DIGIT. 


143 F1C6 


42 


JP $EXACC 




SO GO COLLECT THE EXPONENT. 


144 




$ESIGN: 




145 F1C7 


02 


SET C 




146 F1C8 


6D 


JP JEXCHR 




147 




; 




148 




$EXACC : 




149 




; ACCUMULATE THE EXPLICIT 10 'S EXPONENT. 


150 F1C9 


9100 


LD K, 




151 F1CB 


07 


IF C 
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152 F1CC 


9 IFF 




LD K, OFF 


GET SIGN BITS SET. 




153 F1CE 


AFCA 




PUSH K 


SAVE EXPLICIT EXPONENTS SIGN. 




154 FIDO 


9300 




LD X, 






155 F1D2 


AFCE 




PUSH X 


ZERO EXPLICIT EXPONENT COLLECTED SO 


FAR. 


156 




$EXCLP 








157 F1D4 


02 




SET C 






158 F1D5 


8230C8EB 




SUBC A, '0' 


GET INTEGER FROM ASCII DIGIT. 




159 F1D9 


ACC8CE 




LD X, A 






160 F1DC 


3FC8 




POP A 


GET EXPLICIT EXPONENT COLLECTED SO 


FAR. 


161 FIDE 


A0C8CEF8 




ADD X, A 


X CONTAINS DIGIT + EXP. 




162 F1E2 


A0C8CEF8 




ADD X, A 


X CONTAINS DIGIT + 2*EXP. 




163 F1E6 


E7 




SHL A 






164 F1E7 


E7 




SHL A 






165 F1E8 


E7 




SHL A 


A CONTAINS 8*EXP. 




166 F1E9 


96CEF8 




ADD A, X 


A CONTAINS DIGIT + 10*EXP. 




167 F1EC 


AFC8 




PUSH A 


SAVE BACK ON STACK. 




168 F1EE 


CO 




LDS A, M(B+) 






169 FIEF 


40 




NOP 


GET NEXT CHAR. 




170 FIFO 


9C00 




IFEQ A, 


IS IT A NULL ? 




171 F1F2 


41 




JP $A10EX 


YES SO ADD EXPLICIT AND IMPLICIT 




172 








10 'S EXPONENT. 




173 








GET HERE MEANS IT IS A DIGIT, 




174 F1F3 


7F 




JP $EXCLP 


SO GO BACK AND ACCUMULATE IT. 




175 




; 








176 




$A10EX 








177 




; DONE 


COLLECTING DIGITS. 


ADD THE EXPLICIT AND IMPLICIT 




178 




; 10 'S 


EXPONENT COLLECTEI 


SO FAR. 




179 F1F4 


3FC8 




POP A 


GET EXPLICIT 10 'S EXPONENT. 




180 F1F6 


3FCA 




POP K 


GET IIS SIGN. 




181 F1F8 


8200CAFC 




IFEQ K, 


IS IT POSITIVE ? 




182 F1FC 


42 




JP JADDEX 


YES SO ADD 'EM. 




183 F1FD 


01 




COMP A 






184 FIFE 


04 




INC A 


CHANGE TO 2'S COM. 




185 




$ADDEX 








186 F1FF 


AGFFFAC4F8 




ADD A, W(SP-06) 


ADD IMPLICIT EXPONENT. 




187 F204 


BD7FFF 




IFGT A, 07FFF 


IS IT NEGATIVE ? 




188 F207 


43 




JP JNEG10 


YES, CHANGE IT. 




189 F208 


9300 




LD X, 


LOAD POST. SIGN IN X. 




190 F20A 


44 




JP $ESAVE 






191 




$NEG10 








192 F20B 


93FF 




LD, X, OFF 


LOAD NEG. SIGN IN X. 




193 F20D 


01 




COMP A 






194 F20E 


04 




INC A 


MAKE IT POSITIVE. 




195 




$ESAVE 








196 F20F 


ACC8CC 




LD B, A 


SAVE 10 'S EXPONENT IN A. 




197 F212 


3FC8 




POP A 


GET HI-INT. 




198 F214 


3FCA 




POP K 


GET LO-INT. 




199 F216 


AFCE 




PUSH X 


SAVE SIGN OF 10 ' S EXPONENT. 




200 F218 


AFCC 




PUSH B 


AND ITS VALUE. 




201 




; 








202 




; NOW CONVERT HI-INT. LO-I 


NT TO A NORMALIZED FLOATING POINT 
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AT OF. MAC 










203 




NUMBER. THE BINARY EXPONENT IS COLLECTED IN B. 


204 










205 F21A 


9000 




IFGT A, 


IF HI-INT IS NOT 0, 


206 F21C 


58 




JP $N0RM2 


NEED TO SHIFT K AND A. 


207 F21D 


8200CAFD 




IFGT K, 


IF HI-INT IS 0, BUT NOT LO-INT, 


208 F221 


4E 




JP $N0RM1 


NEED TO SHIFT ONLY K. 


209 








GET HERE MEANS MANTISSA IS 0. 


210 F222 


00 




CLR A 




211 F223 


ACC8CA 




LD K, A 




212 F226 


02 




SET C 




213 F227 


8208C4EB 




SUB SP, 08 


ADJUST SP. DONE! ! ! 


214 F22B 


3FCC 




POP B 




215 F22D 


3FCE 




POP X 




216 F22F 


3C 




RET 




217 










218 


$N0RM1 : 






219 




HI-INT IS 0, SO WORK WITH LO-INT ONLY. 


220 F230 


AECA 




X A, K 




221 F232 


9210 




LD B, 010 


LOAD 16 INTO EXPONENT COUNTER. 


222 F234 


42 




JP pRLUP 




223 










224 


$N0RM2 : 






225 




HI-INT IS NOT 0, SO NEED TO HANDLE BOTH. 


226 F235 


9220 




LD B, 020 


LOAD 32 INTO LOOP COUNTER. 


227 


$NRLUP : 






228 F237 


E7 




SHL A 




229 F238 


07 




IF C 


DID A 1 COME OUT ? 


230 F239 


4D 




JP $NRDUN 


YES IT IS NORMALIZED NOW. 


231 F23A 


AECA 




X A, K 




232 F23C 


E7 




SHL A 




233 F23D 


07 




IF C 




234 F23E 


96CA08 




SET K.O 




235 F241 


AECA 




X A, K 




236 F243 


AACC 




DECSZ B 




237 F245 


40 




NOP 


SHOULD NEVER BE SKIPPED ! ! 


238 F246 


6F 




JP $NRLUP 




239 


$NRDUN : 






240 F247 


D7 




RRC A 


RESTORE SHIFTED 1. 


241 F248 


ABOO 




ST A, TMP1 


STORE IN W(0) . 


242 F24A 


3FCE 




POP X 


GET 10 'S EXPONENT. 


243 F24C 


3FC8 




POP A 


GET 10 'S EXPONENT SIGN. 


244 F24E 


AEOO 




X A, TMP1 


A IS HI-INT ONCE MORE. 


245 F250 


AFCC 




PUSH B 


SAVE BINARY EXPONENT. 


246 F252 


AFCE 




PUSH X 


SAVE 10 'S EXPONENT. 


247 F254 


AECA 




X A, K 


HI-INT TO K, LO-INT TO A. 


248 F256 


960010 




IF TMP1.0 


IS 10 'S EXPONENT NEGATIVE ? 


249 F259 


58 




JP $DIV10 


YES, GO TO DIVIDE BY 10. 


250 




GET HERE MEANS 10 ' S EXI 


ONENT IS POSITIVE, SO MULTIPLY BY 10. 


251 




ACTUALLY, WHAT IS USED 


IS 


252 






10»N = (0.625* (2*4 


)*N 


253 




SO MULTIPLY BY 0.625 NC 


W AND TAKE CARE OF 2*(4*N) LATER. 
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254 F25A 


A4F26ACCAB 




LD B, W($MTLO) 






255 F25F 


AFCC 




PUSH B 


SAVE LO WORD OF 0.625 ON STACK. 




256 F261 


A4F26CCCAB 




LD B, W($MTHI) 






257 F266 


AFCC 




PUSH B 


SAVE HI WORD OF 0.625 ON STACK. 




258 F268 


57 




JP JJAMIT 


GO TO ROUTINE THAT JAMS 0.625"N 




259 








BY REPEATED MULTIPLICATION INTO HI-INT 


LO-INI. 


260 












261 




DEFINE SOME CONSTANTS. 






262 F269 


40 




. EVEN 


FORCE EVEN ADDRESS. 




263 F26A 


0000 $MTLO: 


.WORD 






264 F26C 


OOAO $MTHI : 


.WORD OAOOO 






265 F26E 


CDCC $DTLO: 


.WORD OCCCD 






266 F270 


CCCC fDTHI: 


.WORD OCCCC 






267 












268 


$DIV10: 








269 




GET HERE MEANS 10'S EXPONENT IS NEGATIVE, SO DIVIDE BY 10. 




270 




ACTUALLY WHAT IS DONE IS 




271 






10*(-N) = ((2 A 3)/(0.8)) A (-N) = ((0.8) A N)*(2 A (-3*N))) 




272 




SO MULTIPLY BY 0.8 NOW 


AND TAKE CARE OF 2 A (-3*N) LATER. 




273 F272 


A4F26ECCAB 




LD B, W($DTLO) 






274 F277 


AFCC 




PUSH B 


SAVE LO WORD OF .8 




275 F279 


A4F270CCAB 




LD B, W($DTHI) 






276 F27E 


AFCC 




PUSH B 


SAVE HI WORD OF .8 




277 












278 


I 


JAMIT : 








279 




JAM IN THE MULTIPLICATION PARI NEEDED TO HANDLE THE 10 'S EXP. 




280 F280 


9200 




LD B, 


B IS USED TO TRACK ANY BINARY POWERS 




281 








THAT COME UP DURING NORMALIZATION. 




282 F282 


8200CEFC 




IFEQ X, 


IS 10' S EXPONENT ? 




283 F286 


57 




JP JJAMDN 


YES, DONE ALREADY. 




284 


1 


JAMLP : 








285 F287 


35C0 




JSR BFMUL 


MULTIPLY USING 32 BIT UNSIGNED. 




286 F289 


AECA 




X A, K 


SWAP HI AND LO WORDS. 




287 F28B 


E7 




SHL A 






288 F28C 


07 




IF C 


IS THERE A CARRY ? 




289 F28D 


4A 




JP $ ISNED 


YES, SO IT IS ALREADY NORMALIZED. 




290 F28E 


A9CC 




INC B 


NEED TO SHIFT LEFT TO NORMALIZE, SO 




291 








INCREASE B BY 1. 




292 F290 


AECA 




X A, K 






293 F292 


E7 




SHL A 






294 F293 


07 




IS C 






295 F294 


96CA08 




SET K.O 






296 F297 


43 




JP $0VR1 






297 


1 


ISNED: 








298 F298 


D7 




RRC A 






299 F299 


AECA 




X A, K 






300 


1 


0VR1: 








301 F29B 


AACE 




DECSZ X 


DONE YET ? 




302 F29D 


76 




JP $ JAMLP 


NO SO DO IT ONCE MORE. 




303 




GET HERE MEANS MULTIPLI 


CATIONS HAVE BEEN DONE. NOW TAKE 




304 




CARE 


OF THE EXPONENTS. 







23 



NATIONAL 


SEMICONDUCTOR 


CORPORATION PAGE: 


21 


HPC CROSS ASSEMBLER, REV 


:C,30 JUL 


86 




ATOF 










AT OF. MAC 










305 




$ JAMDN : 






306 F29E 


3FCE 




POP X 




307 F2A0 


3FCE 




POP X 


GET 0.625 OR 0.8 OFF THE STACK. 


308 F2A2 


3FCE 




POP X 


GET THE 10' S EXPONENT. 


309 F2A4 


AFC8 




PUSH A 


SAVE LO WORD OF FLP NUMBER. 


310 F2A6 


AFC A 




PUSH K 


SAVE HI WORD OF FLP NUMBER. 


311 F2A8 


A6FFFAC4A8 




LD A, W(SP-6) 


GET THE BINARY EXPONENT THAT WAS SAVED. 


312 F2AD 


02 




SET C 




313 F2AE 


96CCEB 




SUBC A, B 


SUBTRACT FROM IT BINARY EXPONENT COLLECTED 


314 








DURING THE JAMMING. 


315 F2B1 


ACC8CC 




LD B, A 


SAVE IT IN B. 


316 F2B4 


A8CE 




LD A, X 


MOVE THE 10 'S EXPONENT TO A. 


317 F2B6 


960010 




IF TMP1.0 


IS THE 10 'S EXPONENT NEGATIVE ? 


318 F2B9 


49 




JP $NAGAS 


YES, SO GOT TO SUBTRACT. 


319 








GET HERE MEANS 10' S EXPONENT IS 


320 








POSITIVE, SO MUL IT BY 4. 


321 F2BA 


E7 




SHL A 


MULTIPLY BY 2. 


322 F2BB 


E7 




SHL A 


MULTIPLY BY 2 AGAIN. 


323 F2BC 


96CCF8 




ADD A, B 


GET THE BINARY EXPONENT IN ALSO. 


324 F2BF 


B8007E 




ADD A, 07E 


AND THE IEEE BIAS. 


325 F2C2 


4C 




JP $EXCPT 


GO CHECK FOR OVER/UNDERFLOW. 


326 




; 






327 




$NAGAS : 






328 




; GET HERE MEANS 10 ' S EXPONENT IS NEGATIVE, SO GOT TO MULTIPLY 


329 




; IT BY 


-3. 




330 F2C3 


E7 




SHL A 


MULTIPLY BY 2. 


331 F2C4 


96CEF8 




ADD A, X 


ADD TO GIVE MULTIPLY BY 3. 


332 F2C7 


01 




COMP A 




333 F2C8 


04 




INC A 


MAKE IT NEGATIVE. 


334 F2C9 


96CCF8 




ADD A, B 


GET IN THE BINARY EXPONENT. 


335 F2CC 


B8007E 




ADD A, 07E 


AND THE IEEE BIAS. 


336 




JEXCPT : 






337 




; CHECK 


FOR OVERFLOW/UNDI 


RFLOW. 


338 F2CF 


ACC4CE 




LD X, SP 


FIRST DO SOME JUGGLING 


339 F2D2 


02 




SET C 


TO BE COMPATIBLE WITH EXCEPTION 


340 F2D3 


820ACEEB 




SUBC X, OA 


HANDLING IN OTHER ROUTINES. 


341 F2D7 


AFCE 




PUSH X 




342 F2D9 


BD7FFF 




IFGT A, 07FFF 


IS BIASED EXPONENT NEGATIVE ? 


343 F2DC 


B4FD3F 




JMPL UNDFL 




344 F2DF 


9C00 




IFEQ A, 


IS IT ? 


345 F2E1 


B4FD3A 




JMPL UNDFL 


YES IT IS STILL UNDERFLOW. 


346 F2E4 


9DFE 




IFGT A, OFE 


IS IT GT THAN 254 ? 


347 F2E6 


B4FD46 




JMPL OVRFL 




348 




; GET HERE MEANS VALID SI 


FLP NUMBER. 


349 F2E9 


3FCE 




POP X 


X POINTS TO MANTISSA SIGN. 


350 F2EB 


E7 




SHL A 




351 F2EC 


E7 




SHL A 




352 F2ED 


E7 




SHL A 




353 F2EE 


E7 




SHL A 




354 F2EF 


E7 




SHL A 




355 F2FD 


E7 




SHL A 
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356 


F2F1 


E7 


357 


F2F2 


E7 


358 


F2F3 


8FFA 


359 


F2F5 


ABOO 


360 


F2F7 


3FCA 


361 


F2F9 


3FC8 


362 


F2FB 


Fl 


363 


F2FC 


A8CA 


364 


F2FE 


Fl 


365 


F2FF 


A800 


366 


F301 


F3 


367 


F302 


3664 


368 


F304 


F2 


369 


F305 


F2 


370 


F306 


AFCE 


371 


F308 


368C 


372 


F30A 


3FC4 


373 


F30C 


3FCC 


374 


F30E 


3FCE 


375 


F310 


3C 


376 






377 







SHL A 




SHL A 


MOVE EXPONENT TO HIGH BYTE 


OR A, W(X) 


GET THE MANTISSA SIGN IN. 


ST A, TMP1 


SAVE IT IN TMP1. 


POP K 


FI-HI TO K. 


POP A 


Fl-LO TO A. 


X A, W(X+) 


SAVE Fl-LO. 


LD A, K 




X A, W(X+) 


SAVE Fl-HI. 


LD A, TMP1 




X A, W(X-) 


SAVE F1-EXP.F1-SIGN, X POI 


JSRL SROUND 


ROUND THE RESULT. 


LD A, W(X-) 


X POINTS TO Fl-HI. 


LD A, W(X-) 


X POINTS TO Fl-LO. 


PUSH X 




JSR FPAK 


PACK IT INTO IEEE FORMAT. 


POP SP 




POP B 




POP X 




RET 





POINTS TO Fl-HI. 



.END 
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31 






.FORM 'FTOA. MAC 


32 






.INCLD FTOA. MAC 


1 






.TITLE FTOA 


2 
3 

4 






.LOCAL 






THIS SUBROUTINE CONVERTS A SINGLE PRECISION, BINARY FLOATING 


5 






POINT NUMBER IN THE IEEE FORMAT TO A DECIMAL FLOATING POINT 


6 






STRING. THE DECIMAL FLOATING POINT STRING IS OBTAINED TO A 


7 
8 
9 






PRECISION OF 9 DECIMAL DIGITS. 






THE ALGORITHM USED IS BASED ON: 


10 






J.T. COONEN, 'AN IMPLEMENTATION GUIDE TO A PROPOSED STANDARD 


11 






FOR FLOATING POINT ARITHMETIC,' IEEE COMPUTER, JAN. 1980, PP 68-79. 


12 








13 






ON INPUT, THE BINARY SP FLP NUMBER IS IN REGS. K AND A. 


14 






B CONTAINS THE ADDRESS OF THE LOCATION WHERE THE DECIMAL FLOATING 


15 






POINT STRING IS TO START. NOTE THAT AT LEAST 17 BYTES ARE NEEDED 


16 






FOR THE STORAGE OF THE STRING. THE LAST BYTE IS ALWAYS NULL. 


17 








18 






ALL REGISTERS ARE PRESERVED BY THIS SUBROUTINE. 


19 








20 




FTOA: 


21 


F311 


AFCE 


PUSH X ; SAVE X ON THE STACK. 


22 


F313 


AFCC 


PUSH B ; SAVE B ON THE STACK. 


23 






CHECK AND SEE IF Fl IS A NAN. 


24 


F315 


3605 


JSR FNACHK 


25 


F317 


07 


IF C 


26 


F318 


B401B4 


JMPL $NAN ; YET IT IS, SO GET OUT. 


27 






CHECK AND SEE IF Fl IS ZERO. 


28 


F31B 


36CA 


JSR FZCHK 


29 


F31D 


07 


IF C 


30 


F31E 


B401C8 


JMPL {ZERO ;YES IT IS, SO GET OUT. 


31 






GET HERE MEANS Fl IS A NON-ZERO, NON-NAN FLP NUMBER. 


32 


F321 


ACC4CE 


LD X, SP 


33 


F324 


8206C4F8 


ADD SP, 06 


ADJUST SP. 


34 


F328 


36C7 


JSR FUNPAK 


UNPACK THE NUMBER. 


35 








X POINTS ONE WORD PAST F1-EXP.F1-SIGN 


36 








ON RETURN. 


37 








38 






COMPUTE THE EXPONENT OF 10 FOR DECIMAL FLP NO. 


39 






THIS IS DONE AS FOLLOWS: 


40 






SUPPOSE Fl = FM * (2*11) 


41 






LET U = M*L0G(2) NOTE: LOG IS TO BASE 10. 


42 






THEN V = INT(U+l-9) 


43 






IS USED AS THE 10' S EXPONENT. 


44 






NOTE: INT REFERS TO INTEGER PART. 


45 








46 


F32A 


D2 


LD A, M(X-) 


X POINTS TO Fl-EXP. 


47 


F32B 


D2 


LD A, M(X-) 


LOAD Fl-EXP. X POINTS TO Fl-SIGN. 


48 


F32C 


B7000000 


LD TMP1, 


FIRST GUESS POSITIVE SIGN FOR EXP. 


49 


F330 


B8FF82 


ADD A, 0FF82 


REMOVE IEEE BIAS FROM Fl-EXP. 
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50 F333 


AFC8 


51 F335 


AFCE 


52 F337 


07 


53 F338 


46 


54 F339 


B700FF00 


55 F33D 


01 


56 F33E 


04 


57 




58 




59 F33F 


BE4D10 


60 




61 F342 


AECE 


62 F344 


960010 


63 F347 


41 


64 F348 


4B 


65 




66 F349 


01 


67 F34A 


AECE 


68 F34C 


01 


69 F34D 


B80001 


70 F350 


AECE 


71 F352 


07 


72 F353 


04 


73 




74 F354 


04 


75 F355 


B8FFF7 


76 F358 


BD7FFF 


77 F35B 


45 


78 F35C 


B7000000 


79 F36D 


4F 


80 




81 F361 


B700FF00 


82 F365 


01 


83 F366 


AECE 


84 F368 


01 


85 F369 


B80001 


86 F36C 


AECE 


87 F36E 


07 


88 F36F 


04 


89 




90 




91 




92 




93 




94 




95 




96 




97 




98 




99 




100 





PUSH A 

PUSH X 

IF C 

JP $ML0G2 

LD TMP1, OFF 

COMP A 

INC A 
$ML0G2 : 
; MULTIPLY M BY LOG (2) . 

MULT A, 04D10 

X A, X 

IF TMP1.0 
JP JCSIGN 
JP JREMV9 



$CSIGN: 



COMP A 
X A, X 
COMP A 
ADD A, 
X A, X 
IF C 
INC A 



01 



$REMV9 : 



INC A 
ADD A, 

IFGT A 
JP $CHNGS 
LD TMP1, 
JP $DIV10 



0FFF7 
07FFF 



SAVE IT ON THE STACK. 

SAVE Fl-SIGN ADDRESS ALSO. 

WAS THERE A CARRY ON THE LAST ADD 

YES, SO 2'S EXP IS POSITIVE. 

2'S EXPONENT IS NEGATIVE. 

MAKE IT POSITIVE. 



L0G(2) IS 0.0100110100010000 TO 16 BITS. 

X CONTAINS INTEGER PART, AND A FRACT. PART. 

SWAP THE TWO. 

WAS THE 2'S EXPONENT NEGATIVE ? 

YES, SO MAKE U NEGATIVE. 

NO, SO GO DO V = U + 1 - 9. 



$CHNGS : 



OFF 



COMP INTEGER PART. 



FRACTION PARI. 



INCREASE FRACTION PART. 

SUBTRACT 9. 

IS IT NEGATIVE ? 

YES, SO CHANGE ITS SIGN. 

REMEMBER POSITIVE SIGN. 



REMEMBER NEGATIVE SIGN. 
MAKE V POSITIVE. 



01 



LD TMP1 
COMP A 
X A, X 
COMP A 
ADD A 
X A, X 
IF C 
INC A 
$DIV10: 

V = INI (U+l-9) HAS BEEN COMPUTED AND IS IN A. 

NOW COMPUTE W = F1/(10 A V). W SHOULD BE AN INTEGER, AND IT IS 

COMPUTED TO A 32 BIT PRECISION. 

THIS COMPUTATION IS DONE AS FOLLOWS: 

IF V > 0, THEN F1/(10 A V) = Fl* (0 .8"V) * (2* (-3V) ) . 

IF V < 0, THEN F1/(10 A V) = Fl* (0 .625*U) * (2* (4V) ) . 
SO FIRST MULTIPLY THE MANTISSA OF Fl V TIMES BY 0.8 (OR 0.625) 
AND THEN ADJUST THE EXPONENT OF Fl. NOTE THAT THE PARTIAL PRODUCTS 
IN MULTIPLYING BY 0.8 (OR 0.625) ARE KEPT NORMALIZED. THIS IS 
ESSENTIAL TO PRESERVE 32 BIT ACCURACY IN THE FINAL RESULT. 
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101 




SINCE THE MANTISSA OF Fl IS NORMALIZED, AND 0.8 (OR 0.625 IS ALSO 


102 




NORMALIZED, EACH PRODUCT NEEDS AT MOST 1 LEFT SHIFT FOR 




103 




RENORMALIZATION. THE SHIFTS ACCUMULATED DURING RENORMALIZATION 


ARE 


104 




TRACKED AND ACCOUNTED FOR IN THE CALCULATION. 




105 F370 


3FCE 




POP X 


X NOW POINTS TO Fl-SIGN. 




106 F372 


AFC8 




PUSH A 


SAVE U ON THE STACK. 




107 F374 


ACC8CC 




LD B, A 


MOVE V TO B ALSO. 




108 F377 


F2 




LD A, W(X-) 


X POINTS TO Fl-HI. 




109 F378 


F2 




LD A, W(X-) 


LOAD Fl-HI. X POINTS TO Fl-LO. 




110 F379 


ACC8CA 




LD K, A 






111 F37C 


F4 




LD A, W(X) 


LOAD Fl-LO. 




112 F370 


960010 




IF TMP1.0 


IS V NEGATIVE? 




113 F380 


57 




JP $MUL10 


YES, SO MULTIPLY V TIMES BY .625. 




114 








GET HERE MEANS MULTIPLY V TIMES BY .8. 




115 F381 


A4F390CEAB 




LD X, W($DILO) 






116 F386 


AFCE 




PUSH X 


LO WORD OF 0.8 TO STACK. 




117 F388 


A4F392CEAB 




LD X, W($DTHI) 






118 F38D 


AFCE 




PUSH X 


HI WORD OF 0.8 TO STACK. 




119 F38F 


56 




JP ftJAMIT 


GO DO MULTIPLICATION. 




120 












121 






.EVEN 


FORCE EVEN ADDRESS. 




122 F390 


CDCC JDTLO 


.WORD OCCCD 






123 F392 


CCCC $DTHI 


.WORD OCCCC 






124 F394 


0000 $MILO 


.WORD 






125 F396 


OOAO $MIHI 


.WORD OAOOO 






126 












127 


$MUL1C 


: 






128 F398 


A4F394CEAB 




LD X, W($MILO) 






129 F39D 


AFCE 




PUSH X 


LO WORD OF 0.625 TO STACK. 




130 F39F 


A4F396CEAB 




LD X, W($MTHI) 






131 F3A4 


AFCE 




PUSH X 


HI WORD OF 0.625 TO STACK. 




132 












133 


$JAMI1 








134 F3A6 


9300 




LD X, 


INIT X TO TRACK ANY POWERS OF 




135 








2 GENERATED DURING NORMALIZATION 




136 








OF PARTIAL PRODUCTS. 




137 F3A8 


8200CCFC 




IFEQ B, 


IS B ALREADY ? 




138 F3AC 


57 




JP ftJAMON 


YES, SO SKIP MULTIPLY LOOP. 




139 


$JAMLI 


3 . 






140 F3AD 


36E6 




JSR BFMUL 


MULTIPLY. 




141 F3AF 


AECA 




X A, K 


SWAP HI AND LO WORDS OF PART. PROD. 




142 F3B1 


E7 




SHL A 






143 F3B2 


07 




IF C 


IS THERE A CARRY ? 




144 F3B3 


4A 




JP $ISNED 


YES, SO SKIP OVER RENORMALIZATION. 




145 








GET HERE MEANS NEED TO RENORM. 




146 F3B4 


A9CE 




INC X 


UPDATE RENORM COUNT. 




147 F3B6 


AECA 




X A, K 


SWAP HI AND LO PART. PROD. 




148 F3B8 


E7 




SHL A 






149 F3B9 


07 




IF C 






150 F3BA 


96CA08 




SET K.O 


SET BIT SHIFTED OUT FROM LO WORD. 




151 F3BD 


43 




JP $0VR1 
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152 




$ISNED: 








153 F3BE 


D7 




RRC A 




PUT BACK SHIFTED BIT. 


154 F3BF 


AECA 




X A, K 






155 




$0VR1 : 








156 F3C1 


AACC 




DECSZ B 




IS B YET ? 


157 F3C3 


76 




JP $JAMLP 




NO, SO DO IT AGAIN. 


158 




$JAMON: 








159 




; GET HERE MEANS MULTIPLICATION HAS BEEN DONE, SO TAKE CARE 


160 




; OF EXPONENT. 






161 F3C4 


3FCC 




POP B 






162 F3C6 


3FCC 




POP B 




GET RID OF 0.8 (OR 0.625) FROM STACK. 


163 F3C8 


AFC8 




PUSH A 




SAVE LO WORD OF PROD. 


164 F3CA 


AFC A 




PUSH K 




SAVE HI WORD OF PRODUCT. 


165 F3CC 


A6FFF8C4A8 




LD A, W(SP- 


38) 


GET Fl'S BINARY EXPONENT. 


166 F3D1 


02 




SET C 






167 F3D2 


96CEEB 




SUBC A, X 




SUBTRACT FROM IT ANYTHING COLLECTED 


168 










DURING RENORM. 


169 F3D5 


ACC8CE 




LD X, A 




AND SAVE IT IN X. 


170 F3D8 


A6FFFAC4A8 




LD A, W(SP- 


36) 


GET V FROM THE STACK. 


171 F3DD 


960010 




IF TMP1.0 




IS V NEGATIVE ? 


172 F3E0 


49 




JP $ML4 




YES, SO MULTIPLY V BY 4. 


173 










GET HERE MEANS MULTIPLY V BY -3. 


174 F3E1 


E7 




SHL A 




NOW A CONTAINS 2*V. 


175 F3E2 


A6FFFAC4F8 




ADD A, W(SP 


-06) 


NOW A CONTAINS 3*V. 


176 F3E7 


01 




COMP A 






177 F3E8 


04 




INC A 




NOW A CONTAINS -3*V. 


178 F3E9 


42 




JP JADEM 




GO FIGURE FINAL EXPONENT. 


179 




$ML4: 








180 F3EA 


E7 




SHL A 






181 F3EB 


E7 




SHL A 




NOW A CONTAINS 4*V. 


182 




$ADEM : 








183 F3EC 


96CEF8 




ADD A, X 




A SHOULD NOW BE A POSITIVE INTEGER 


184 










IN THE RANGE TO 32. 


185 




; NOW CHECK AND SEE 


IF A 


HAS ENOUGH PRECISION. 


186 F3EF 


9020 




IFGT A, 020 




NEED MORE THAN 32 BITS ? 


187 F3F1 


5A 




JP $INCRV 




YES, SO GO INCREASE V. 


188 F3F2 


9D1B 




IFGT A, 01B 




NEED AT LEAST 28 BITS ? 


189 F3F4 


9435 




JMP $G00N 




YES, SO ALL IS OK. GO ON. 


190 










GET HERE MEANS NEED MORE 


191 










PRECISION, SO DECREASE V. 


192 F3F6 


3FC8 




POP A 




GET HI-PROD OFF STACK. 


193 F3F8 


3FC8 




POP A 




GET LO WORD OFF STACK. 


194 F3FA 


3FC8 




POP A 




GET MAGN. OF V. 


195 F3FC 


96D010 




IF TMP1.0 




IS V NEG. ? 


196 F3FF 


56 




JP $VUP 




YES, SO GO INCR. MAGN. OF V. 


197 










GET HERE MEANS V IS POSITIVE, 


198 










AND NEED TO DECREMENT IT. 


199 F400 


B8FFFF 




ADD A, OFFFF 


SUBTRACT 1 FROM A. 


200 F403 


07 




IF C 




GOT A CARRY ? 


201 F404 


5A 




JP $GOBAK 




THEN OK. 


202 F405 


01 




COMP A 
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203 


F406 


04 


204 


F407 


B700FF00 


205 


F40B 


53 


206 






207 


F40C 


3FC8 


208 


F40E 


3FC8 


209 


F410 


3FC8 


210 


F412 


960010 


211 


F415 


42 


212 






213 


F416 


04 


214 


F417 


47 


215 






216 


F418 


AAC8 


217 


F41A 


44 


218 


F41B 


B7000000 


219 






220 


F41F 


ACC4CE 


221 


F422 


02 


222 


F423 


8204CEEB 


223 


F427 


AFCE 


224 


F429 


95B9 


225 






226 


F42B 


01 


227 


F42C 


04 


228 


F42D 


B80020 


229 


F430 


ACC8CE 


230 


F433 


3FCA 


231 


F435 


3FC8 


232 


F437 


8200CEFC 


233 


F43B 


49 


234 






235 






236 


F43C 


AECA 


237 


F43E 


C7 


238 


F43F 


AECA 


239 


F441 


D7 


240 


F442 


AACE 


241 


F444 


68 


242 






243 






244 






245 


F445 


AFC8 


246 


F447 


AFC A 


247 


F449 


A6FFF0C4A8 


248 


F44E 


B80010 


249 


F451 


ACC8CC 


250 


F454 


00 


251 


F455 


C3 


252 


F456 


40 


253 


F457 


A6FFFAC4A8 



$INCRV: 



$VUP: 



$VDOWN : 



$GOBAK: 



INC A 


LD TMP1, OFF 


JP |GOBAK 


POP A 


POP A 


POP A 


IF TMP1.0 


JP ftVDOWN 


INC A 


JP ftGOBAK 


DECSZ A 


JP $GOBAK 


LD TMP1, 


LD X, SP 


SET C 


SUBC X, 04 


PUSH X 


JMP $DIV10 


COMP A 


INC A 


ADD A, 020 


LD X, A 


POP X 


POP A 


IFEQ X, 


JP $INDUN 


JUST THE PRO 


X A, K 


SHR A 


X A, K 


RRC A 


DECSZ X 


JP $INTFY 



U CHANGES SIGN. 



GET HI PROD. OFF STACK. 
GET LO PROD. OFF STACK. 
GET MAGN. OF V. 
IS V NEGATIVE ? 
YES. 



V CHANGES SIGN. 



$GOON: 

COMP A 

NEGATE A. 

SUBTRACT IT FROM 32. 

AND MOVE IT TO X. 

GET HI WORD OF PRODUCT. 

GET LO WORD OF PROD. 

IS X ? 

YES, SO ALREADY A 32 BIT INTEGER. 
$INTFY: 

TO FORM A 32 BIT INTEGER. 

SWAP HI AND LO WORDS. 



SHIFT IT RIGHT ONCE. 

X YET ? 

NO SO GO DO SOME MORE. 
$INDUN: 

; GET HERE MEANS K.A CONTAIN THE 32 BIT INTEGER THAT IS THE 
; MANTISSA OF THE DECIMAL FLP NUMBER. 



PUSH A 




SAVE LO-INT. 


PUSH K 




SAVE HI INT. 


LD A, W(SP- 


-010) 


GET STARTING ADDRESS OF DECIMAL STRING 


ADD A, 010 




ADD 16 TO IT. 


LD B, A 




AND MOVE IT B. 


CLR A 






XS A, M(B- 




OUTPUT TERMINATING NULL BYTE. 


NOP 






LD A, W(SP- 


■06) 


GET V. 
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254 


F45C 


9F0A 


255 






256 


F45E 


AECE 


257 


F460 


B80030 


258 


F463 


C3 


259 


F464 


40 


260 


F465 


A8CE 


261 


F467 


B80030 


262 


F46A 


C3 


263 


F46B 


40 


264 


F46C 


902B 


265 


F46E 


960010 


266 


F471 


902D 


267 


F473 


C3 


268 


F474 


40 


269 


F475 


9045 


270 


F477 


C3 


271 


F478 


40 


272 


F479 


902E 


273 


F47B 


C3 


274 


F47C 


40 


275 






276 


F47D 


B7000A00 


277 






278 


F481 


3FC8 


279 


F483 


9F0A 


280 






281 


F485 


ACC8CA 


282 


F488 


3FC8 


283 


F48A 


AFCC 


284 


F48C 


ACCACC 


285 


F48F 


82 




F490 


OA 




F491 


C8 




F492 


EF 


286 






287 






288 






289 






290 


F493 


ACCCCA 


291 


F496 


3FCC 


292 


F498 


AFC8 


293 


F49A 


AFC A 


294 


F49C 


A8CE 


295 


F49E 


B80030 


296 


F4A1 


C3 


297 


F4A2 


40 


298 


F4A3 


AAOO 


299 


F4A5 


9524 


300 






301 


F4A7 


3FC8 



DIV A, OA 

X A, X 
ADD A, 030 

XS A, M(B-) 

NOP 

LD A, X 

ADD A, 030 

XS A, M(B-) 

NOP 

LD A, 028 

IF TMP1.0 

LD A, 02D 

XS A, M(B-) 

NOP 

LD A, 045 

XS A, M(B-) 

NOP 

LD A, 02E 

XS A, M(B-) 

NOP 

; NOW NEED TO OUTPUT 
LD TMP1, OA 

$DOLUP: 

POP A 
DIV A, OA 

LD K, A 
POP A 
PUSH B 
LD B, K 
.BYTE 082, OA. 



DIVIDE IT BY 10. QUOT. IN A, 

REM. IN X. 

REM TO A. 

MAKE IT INTO ASCII BYTE. 

OUTPUT IT. 



FINISHED OUTPUTING EXPONENT. 
SAY EXP SIGN IS '+' . 



; NOPE, IT IS '-' 
; OUTPUT IT. 



OUTPUT 'E' 



OUTPUT 



10 DECIMAL DIGITS. 

: LOAD 10 INTO TMP1 AS LOOP COUNTER. 



A CONTAINS HI INT. 

DIVIDE IT BY 10. QUOT. IN A, 

REM. IN X. 



A CONTAINS LO INT. 
SAVE DEC. STR. ADDR. 
B CONTAINS HI-QUOT. 
0C8,0EF 



THE ABOVE 4 BYTES REPRESENT THE INSTRUCTION DIVD A, OA. 
BECAUSE THE ASSEMBLER DOES NOT KNOW ABOUT IT YET, WE HAVE TO 
KLUDGE IT THIS WAY. 
AFTER THE DIVD, A CONTAINS THE LO-QUOT. AND X THE REM. 



LD K, B 
POP B 
PUSH A 
PUSH K 
LD A, X 
ADD A, 030 
XS A, M(B-) 
NOP 

DECSZ TMP1 
JMP $DOLUP 



MOVE HI-QUOT TO 
B CONTAINS DEC. 
SAVE LO INI. 
SAVE HI INT. 
MOVE REM TO A. 
ASCII-FY II. 
AND OUTPUT IT. 



AND 

K. 

STR. 



ADDR. 



IS 

NO. 



TMP1 YET ? 
GO GET SOME MORE. 



GET 



HERE MEANS 
POP A 



DONE WITH OUTPUTING MANTISSA. 
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302 F4A9 


3FC8 




POP A 




303 F4AB 


3FC8 




POP A 




304 F4AD 


3FC8 




POP A 


GET SOME GARBAGE OFF THE STACK. 


305 F4AF 


3FCA 




POP K 


GET F1-EXP.F1-SIGN TO K. 


306 F4B1 


9020 




LD A, 02D 


LOAD SP INTO A. 


307 F4B3 


96CA10 




IF K.O 




308 F4B6 


902D 




LD A, 02D 


IF MAINT. IS NEG. LOAD '-'. 


309 F4B8 


C6 




ST A, M(B) 


OUTPUT SIGN. 


310 F4B9 


AFC A 




PUSH K 


F1-EXP.F1-SIGN BACK ON STACK. 


311 F4BB 


ACC4CE 




LD X, SP 




312 F4BE 


02 




SET C 




313 F4BF 


8206CEEB 




SUBC X, 06 


X POINTS TO Fl-LO. 


314 F4C3 


AFCE 




PUSH X 




315 F4C5 


B5FBB4 + 




JSR FPAK 


PACK IT, SO RESTORING K AND A. 


316 F4C8 


3FC4 




POP SP 




317 F4CA 


3FCC 




POP B 


RESTORE B. 


318 F4CC 


3FCE 




POP X 


RESTORE X. 


319 F4CE 


3C 




RET 




320 




; 






321 




$NAN: 






322 




; GET HERE MEANS Fl IS A 


NAN. 


323 F4CF 


AFC8 




PUSH A 




324 F4D1 


ACCCCE 




LD X, B 




325 F4D4 


8210CEF8 




ADD X, 010 




326 F4D8 


00 




CLR A 




327 F4D9 


03 




X A, M(X-) 




328 F4DA 


9210 




LD B, 010 




329 




JNANLP : 






330 F4DC 


9 OFF 




LD A, OFF 




331 F4DE 


D3 




X A, M(X-) 




332 F4DF 


AACC 




DECSZ B 




333 F4E1 


65 




JP $NANLP 




334 F4E2 


3FC8 




POP A 




335 F4E4 


3FCC 




POP B 




336 F4E6 


3FCE 




POP X 




337 F4E8 


3C 




RET 




338 




; 






339 




$ZERO: 






340 




; GET HERE MEANS Fl IS ZERO. 


341 F4E9 


AFC8 




PUSH A 




342 F4EB 


ACCCCE 




LD X, B 


X CONTAINS DECIMAL STRING ADDR. 


343 F4EE 


8210CEF8 




ADD X, 010 




344 F4F2 


00 




CLR A 




345 F4F3 


D3 




X A, M(X-) 


OUTPUT TERMINATING NULL BYTE. 


346 F4F4 


9030 




LD A, 030 


LOAD INTO A. 


347 F4F6 


D3 




X A, M(X-) 




348 F4F7 


9030 




LD A, 030 




349 F4F9 


D3 




X A, M(X-) 


OUTPUT 00 FOR EXPONENT. 


350 F4FA 


902B 




LD A, 02B 


LOAD '+' SIGN. 


351 F4FC 


D3 




X A, M(X-) 




352 F4FD 


9045 




LD A, 045 


LOAD 'E' 
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353 F4FF D3 

354 F500 902E LD A, 02E : LOAD 

355 F502 D3 

356 F503 920A 

357 $ZERLP: 

358 F505 9030 

359 F507 D3 

360 F508 AACC 

361 F50A 65 

362 F508 9020 LD A, 020 : LOAD SP. 

363 F50D D5 

364 F50E 3FC8 

365 F510 3FCC 

366 F512 3FCE 

367 F514 3C 

368 ; 

369 .END 



X A, 


M(X-) 


LD A 


02E 


X A, 


M(X-) 


LD B 


OA 


LD A 


030 


X A, 


M(X-) 


DECSZ B 


JP JZERLP 


LD A 


020 


X A, 


M(X) 


POP A 


POP B 


POP X 


RET 
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33 






.FORM, 'FADD. MAX' 


34 






.INCLD FADD. MAC 


1 






.TITLE FADD 


2 
3 

4 






.LOCAL 




SUBROUTINE TO ADD/SUBTRACT TWO SP FLOATING POINT NUMBERS. 


5 
6 
7 






C = Fl + F2 OR C = Fl - F2 




Fl IS 


STORED IN THE IEEE FORMAT IN REGS K AND A. 


8 




THE HIGH WORD OF Fl WILL BE REFERRED AS Fl-Rl AND IS IN K. 


9 




THE LOW WORD OF Fl WILL BE REFERRED TO AS Fl-RO AND IS IN A. 


10 








11 




F2 IS 


STORED IN THE IEEE FORMAT ON THE STACK. IF SP IS THE 


12 




STACK 


POINTER ON ENTRY, THEN 


13 




THE HIGH WORD OF F2 , REFERRED TO AS F2-R1 IS AT SP - 4 AND 


14 




THE LOW WORD OF F2, REFERRED TO AS F2-R0 IS AT SP - 6. 


15 








16 




C IS 


RETURNED IN THE IEEE FORMAT IN REGS K AND A. 


17 








18 


I 


SUB: 




19 


F515 ABOO 




ST A, TMP1 


20 


F517 A6FFFAC4A8 




LD A, W(SP-06) 


LOAD F2-R0. 


21 


F51C AFC8 




PUSH A 


AND SAVE ON STACK. 


22 


F51E A6FFFAC4A8 




LD A, W(SP-06) 


LOAD F2-R1. 


23 


F523 BB8000 




XOR A, 08000 


CHANGE THE SIGN. 


24 


F526 AFC8 




PUSH A 


AND SAVE ON THE STACK. 


25 


F528 A800 




LD A, TMP1 


RESTORE A. 


26 


F52A 3009 




JSR FADD 


CALL THE ADD ROUTINE. 


27 


F52C ABOO 




ST A, TMP1 


SAVE A. 


28 


F52E 3FC8 




POP A 


GET RID OF JUNK 


29 


F530 3FC8 




POP A 


FROM THE STACK. 


30 


F532 A800 




LD A, TMP1 


RESTORE A. 


31 


F534 3C 




RET 


32 








33 


I 


ADD: 




34 




SAVE 


ADDRESS OF F2-R0 IN IMP1. 


35 


F535 AFCE 




PUSH X ; SAVE X ON ENTRY. 


36 


F537 AFCC 




PUSH B ; AND B ON ENTRY. 


37 


F539 ACC4CE 




LD X, SP 


38 


F53C 86FFF6CEF8 




ADD X, 0FFF6 ; SUBTRACT 10. 


39 


F541 ACCEOO 




LD TMP1, X ; AND SAVE IN TMP1. 


40 




CHECK 


AND SEE IF Fl IS A NAN. 


41 


F544 B5FAF9 + 




JSR FNACHK 


42 


F547 07 




IF C 


43 


F548 B4FAC4 




JMPL FNAN ; Fl IS A NAN. 


44 




CHECK 


AND SEE IF F2 IS A NAN. 


45 


F54B ACCACC 




LD B, K 


46 


F54E ACC8CE 




LD X, A 


47 


F551 A20200A8 




LD A, W(TMPl+2) 


48 


F555 ACC8CA 




LD K, A 


49 


F558 AECE 




X A, X 
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50 


F55A 


B5FAE3 


51 


F55D 


07 


52 


F55E 


B4FAAE 


53 






54 


F561 


B5FAED 


55 


F564 


06 


56 


F565 


48 


57 


F566 


ACCCCA 


58 


F569 


3FCC 


59 


F56B 


3FCE 


60 


F56D 


3C 


61 






62 






63 


F56E 


ACCCCA 


64 


F571 


B5FADD 


65 


F574 


06 


66 


F575 


4F 


67 


F567 


A20200AB 


68 


F57A 


ACC8CA 


69 


F57D 


AD00A8 


70 


F580 


3FCC 


71 


F582 


3FCE 


72 


F584 


3C 


73 






74 






75 






76 


F585 


ACC4CE 


77 


F588 


8210C4F8 


78 


F58C 


AFCE 


79 


F58E 


B5FAD0 


80 


F591 


ACOOCC 


81 


F594 


ACCEOO 


82 


F597 


EO 


83 


F598 


40 


84 


F599 


AECA 


85 


F59B 


E4 


86 


F59C 


AECA 


87 


F59E 


B5FAC0 


88 






89 


F5A1 


F2 


90 


F5A2 


ACOOCC 


91 


F5A5 


E2 


92 


F5A6 


40 


93 






94 


F5A7 


F2 


95 


F5A8 


B9FF00 


96 


F5AB 


A6FFFCC4AB 


97 


F5B0 


E2 


98 


F5B1 


40 


99 


F5B2 


B9FF00 


100 


F5B5 


02 



+ JSR FNACHK 

IF C 
JMPL FNAN 

; CHECK AND SEE IF F2 I 
+ JSR FZCHK 

IFN C 
JP $F1CHK 
LD K, B 
POP B 
POP X 
RET 

; CHECK AND SEE IF Fl I 
$F1CHK: 

LD K, B 
+ JSR FZCHK 

IFN C 

JP $NTZERO 
LD A, W(TMPl+2) 
LD K, A 
LD A, W(TMP1) 
POP B 
POP X 
RET 

; GET HERE MEANS NORMAL 

; UNPACK Fl AND F2. 
$NTZERO: 

LD X, SP 
ADD SP, 010 
PUSH X 
+ JSR FUNPAK 

LD B, TMP1 
LD TMP1, X 
LDS A, W(B+) 
NOP 

X A, K 
LD A, W(B) 
X A, K 
+ JSR FUNPAK 

; SET X TO POINT TO F2> 
LD A, W(X-) 
LD B, TMP1 
LDS A, W(B-) 
NOP 

; COMPARE Fl-EXP AND F2 
LD A, W(X-) 
AND A, OFFOO 
ST A, W(SP-4) 
LDS A, W(B-) 
NOP 

AND A, OFFOO 
SET C 



; F2 IS NAN. 
S ZERO. 



F2 IS NOT ZERO. CHECK Fl. 
F2 IS ZERO, SO ANSWER IS Fl. 



S ZERO. 
; RESTORE Fl-Rl FROM B. 



JUMP SINCE Fl IS ALSO NOT ZERO. 
GET HERE MEANS Fl IS ZERO, 
SO ANSWER IS F2. 



ADDITION. 



X POINTS TO Fl-LO. 

MOVE SP PAST LOCAL STORAGE. 

SAVE SP ON STACK FOR QUICK RETURN. 

UNPACK Fl. 

B NOW POINTS TO F2-R0 . 

TMP1 NOW POINTS TO F2-L0. 

LOAD F2-R0 INTO A. 



; LOAD F2-R1 INTO K. 
; UNPACK F2. 
SIGN AND B TO POINT TO Fl-SIGN. 



EXP. 

LOAD F2-EXP.F2-SIGN INTO A. 

MASK OUT SIGN. 

SAVE IN C-SIGN.C-EXP. 

LOAD Fl-EXP. Fl-SIGN INTO A. 
CHANGE TO Fl-EXP. 00000000 . 
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101 F5B6 


A6FFFCC4EB 


102 F5BB 


06 


103 F5BC 


942D 


104 




105 F5BE 


80C9CAAB 


106 F5C2 


A6FFFCC4FB 


107 F5C7 


A6FFFCC4AB 


108 F5CC 


8217CAFD 


109 F5D0 


51 


110 




111 F5D1 


8200CAFC 


112 F5D5 


943B 


113 




114 F5D7 


F4 


115 F5D8 


C7 


116 F5D9 


F3 


117 F5DA 


F4 


118 F5DB 


D7 


119 F5DC 


Fl 


120 F5DD 


AACA 


121 F5DF 


68 


122 F5E0 


9430 


123 




124 




125 F5E2 


FO 


126 F5E3 


00 


127 F5E4 


F3 


128 F5E5 


00 


129 F5E6 


F3 


130 F5E7 


00 


131 F5E8 


Fl 


132 F5E9 


9427 


133 




134 




135 F5EB 


01 


136 F5EC 


04 


137 F5ED 


80C9CAAB 


138 F5F1 


8217CAFD 


139 F5F5 


51 


140 




141 F5F6 


8200CAFC 


142 F5FA 


57 


143 




144 F5FB 


E4 


145 F5FC 


C7 


146 F5FD 


E3 


147 F5FE 


40 


148 F5FF 


E4 


149 F600 


D7 


150 F601 


El 


151 F602 


40 



SUBC A, W(SP-4) ; SUBTRACT F2-EXP. 00000000. 

IFN C 

JMP $F2GTR ; F2-EXP IS BIGGER THAN Fl-EXP. 
; GET HERE MEANS Fl-EXP IS BIGGER THAN F2-EXP. 

LD K, H(A) ; SAVE DIFF. IN K TO BE USED AS LOOP COUNTER. 

ADD A, W(SP-4) 

ST A, W(SP-4) ; RESTORE Fl-EXP AND STORE IN C-SIGN. 

IFGT K, 017 

JP $ZR0F2 ; K GT 23-DEC MEANS F2 GETS ZEROED IN SHIFTS. 
; LOOP TO SHIFT F2 INTO ALIGNMENT. 

IFEQ K, 

; K = MEANS DONE SHIFTING. 



JMP $ADDMN 



$L00P2 : 



X POINTS TO F2-EXP.F2-SIGN. 
AND STORE IT BACK. 



LD A, W(X) 

SHR A 

X A, W(X-) 

LD A, W(X) 

RRC A 

X A, W(X+) 

DECSZ K 

JP $L00P2 

JMP $ADDMN 
$ZR0F2: 
; SET F2 MANTISSA TO 0. 

LD A, W(X+) 

CLR A 

X A, W(X-) 

CLR A 

X A, W(X-) 

CLR A 

X A, W(X+) 

JMP $ADDMN 
; F2 EXPONENT IS GREATER THAN Fl EXPONENT. 
$F26TR: 

COMP A 

INC A ; CHANGE DIFF IN EXP TO POSITIVE. 

LD K, H(A) ; LOAD K WITH LOOP COUNTER. 

IFGT K, 017 

JP $ZR0F1 ; Fl MANT. REDUCED TO IN SHIFTS. 
; LOOP TO SHIFT Fl MANT INTO ALIGNMENT. 

IFEQ K, 

JP $ADDMN ; K=0 MEANS DONE SHIFTING. 



$L00P1: 



LD A, W(B) 

SHR A 

XS A, W(B-) 

NOP 

LD A, W(B) 

RRC A 

XS A, W(B+) 

NOP 
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152 F603 


AACA 




DECSZ K 




153 F605 


6A 




JP JL00P1 




154 F606 


4B 




JP $ADDMN 




155 




fZROFl: 




SET Fl MANT TO 0. 


156 F607 


EO 




LOS A, W(B+) 


B POINTS TO F1-EXP.F1-SIGN. 


157 F608 


40 




NOP 




158 F609 


00 




CLR A 




159 F60A 


E3 




XS A, W(B-) 


STORE IT BACK. 


160 F60B 


40 




NOP 




161 F60C 


00 




CLR A 




162 F60D 


E3 




XS A, W(B-) 




163 F60E 


40 




NOP 




164 F60F 


00 




CLR A 




165 F610 


El 




XS A, W(B+) 




166 F611 


40 




NOP 




167 




; DETERMINE IF MANTISSAS 


ARE TO BE ADDED OR SUBTRACTED. 


168 




$ADDMN : 




B POINTS TO Fl-HI, X TO F2-HI. 


169 F612 


EO 




LDS A, W(B+) 




170 F613 


40 




NOP 




171 F614 


FO 




LD A, W(X+) 




172 F615 


D4 




LD A, M(X) 


LOAD F2-SIGN. 


173 F616 


D8 




XOR A, M(B) 


XOR WITH Fl-SIGN. 


174 F617 


9C00 




IFEQ A, 




175 F619 


9451 




JMP $TRADD 


SAME SIGN SO GO TO ADD MANTISSA. 


176 




; GET HERE MEANS TRUE SUBTRACT OF MANTISSA. 


177 F61B 


F2 




LD A, W(X-) 




178 F61C 


F2 




LD A, W(X-) 


X POINTS TO F2-L0. 


179 F61D 


E2 




LDS A, W(B-) 




180 F61E 


40 




NOP 




181 F61F 


E2 




LDS A, W(B-) 




182 F620 


40 




NOP 


B NOW POINTS TO Fl-LO. 


183 F621 


EO 




LDS A, W(B+) 




184 F622 


40 




NOP 


A NOW CONTAINS Fl-LO. 


185 F623 


02 




SET C 




186 F624 


8FEB 




SUBC A, W(X) 


SUBTRACT F2-L0. 


187 F626 


Fl 




X A, W(X+) 




188 F627 


EO 




LDS A, W(B+) 


A CONTAINS Fl-HI. 


189 F628 


40 




NOP 




190 F629 


8FEB 




SUBC A, W(X) 


SUBTRACT F2-HI. 


191 F62B 


Fl 




X A, W(X+) 




192 F62C 


07 




IF C 




193 F62D 


55 




JP $F1SIN 


Fl GE F2, SO SIGN IS Fl-SIGN. 


194 




; GET HERE MEANS Fl LT F2 


, SO SIGN IS F2-SIGN. 


195 F62E 


A6FFFCC4A8 




LD A, W(SP-4) 




196 F633 


8FDA 




OR A, M(X) 




197 F635 


F3 




X A, W(X-) 


C-EXP.C-SIGN HAS BEEN DETERMINED. 


198 F636 


F4 




LD A, W(X) 




199 F637 


Dl 




COMP A 




200 F638 


F3 




X A, W(X-) 




201 F639 


F4 




LD A, W(X) 




202 F63A 


01 




COMP A 
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203 F63B 


B80001 




ADD A, 01 






204 F63E 


Fl 




X A, W(X+) 






205 F63F 


07 




IF C 






206 F640 


8FA9 




INC W(X) 






207 F642 


47 




JP JANORM 






208 




; GET HERE MEANS Fl GE FS 






209 




$F1SIN: 








210 F643 


A6FFFCC4A8 




LD A, W(SP-4) 






211 F648 


DA 




OR A, M(B) 






212 F649 


F3 




X A, W(X-) 






213 




; NORMALIZE THE MANTISSA. 






214 




JANORM : 








215 F64A 


ACCECC 




LD B, X 


B POINTS TO C-HI. 




216 F64D 


EO 




LDS A, W(B+) 






217 F64E 


40 




NOP 






218 F64F 


CO 




LDS A, M(B+) 






219 F650 


40 




NOP 


B NOW POINTS TO C-EXP BYTE. 




220 F651 


9118 




LD K, 018 


SET UP LOOP LIMIT OF 24-DEC 


IN K. 


221 




$NLOOP : 








222 F653 


F4 




LD A, W(X) 






223 F654 


E7 




SHL A 






224 F655 


07 




IF C 


CARRY MEANS NORMALIZED. 




225 F656 


9448 




JMP {ROUND 


SO JUMP TO ROUNDING CODE. 




226 F658 


F3 




X A, W(X-) 






227 F659 


F4 




LD A, W(X) 






228 F65A 


E7 




SHL A 






229 F65B 


Fl 




X A, W(X+) 






230 F65C 


07 




IF C 






231 F65D 


8F08 




SET W(X) .0 






232 F65F 


ADCC8A 




DECSZ M(B) 


ADJUST EXPONENT. 




233 F662 


43 




JP $0V1 






234 F663 


B4F9B8 




JMPL UNDFL 


C-EXP ZERO MEANS UNDERFLOW. 




235 




$0V1: 








236 F666 


AACA 




DECSZ K 


DECREMENT LOOP COUNTER. 




237 F668 


75 




JP JNLOOP 


GO BACK TO LOOP. 




238 F669 


B4F9B2 




JMPL UNDFL 


UNDERFLOW 




239 




;GET HERE MEANS TRUE ADDI 


HON OF MANTISSA. 




240 




$TRADD : 








241 F66C 


E2 




LDS A, W(B-) 






242 F66D 


40 




NOP 






243 F66E 


E2 




LDS A, W(B-) 






244 F66F 


40 




NOP 


B NOW POINTS TO Fl-HI. 




245 F670 


F2 




LD A, W(X-) 






246 F671 


F2 




LD A, W(X-) 






247 F672 


F4 




LD A, W(X) 


LOAD F2-L0 INTO A. 




248 F673 


F8 




ADD A, W(B) 


ADD Fl-LO. 




249 F674 


Fl 




X A, W(X+) 


STORE IN F2-L0. 




250 F675 


EO 




LDS A, W(B+) 






251 F676 


40 




NOP 


B NOW POINTS TO Fl-HI. 




252 F677 


F4 




LD A, W(X) 


LOAD F2-HI INTO A. 




253 F678 


E8 




ADC A, W(B) 


ADD Fl-HI WITH CARRY FROM LO 


ADD. 
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254 F679 


07 


IF C 




255 F67A 


4A 


JP JADJEX 


IF CARRY, NEED TO INCREASE EXP. 


256 F67B 


Fl 


X A, W(X+) 


STORE RESULT IN F2-HI . 


257 F67C 


A6FFFCC4A8 


LD A, W(SP-4) 


GET C-EXP. 00000000. 


258 F681 


8FDA 


OR A, M(X) 


INTRODUCE SIGN. 


259 F683 


F3 


K A, W(X-) 


STORE IN F2-EXP.F2-SIGN. 


260 F684 


5B 


JP $R0UND 




261 




; GET HERE MEANS NEED TO 


INCREASE EXP BY 1. 


262 




$AD JEX : 




263 F685 


D7 


RRC A 




264 F686 


F3 


X A, W(X-) 




265 F687 


F4 


LD A, W(X) 




266 F688 


D7 


RRC A 




267 F689 


Fl 


X A, W(X+) 




268 F68A 


FO 


LD A, W(X+) 


X NOW POINTS TO F2-EXP.F2-SIGN. 


269 F68B 


A6FFFCC4A8 


LD A, W(SP-4) 


GET C-EXP. 00000000. 


270 F690 


B80100 


ADD A, 0100 


INCREASE EXP BY 1. 


271 F693 


07 


IF C 




272 F694 


B4F998 


JMPL OVRFL 




273 F697 


BDFEFF 


IFGT A, OFEFF 


IS BIASED EXPONENT 255-DEC ? 


274 F69A 


B4F992 


JMPL OVRFL 




275 F69D 


8FDA 


OR A, M(X) 




276 F69F 


F3 


X A, W(X-) 




277 




; NEED TO ROUND THE RESULT. X POINTS TO C-HI. 


278 




$ROUND : 




279 F6A0 


B5F9FB 


JSRL SROUND 




280 




; FINAL CHECK OF EXPONENI 




281 F6A3 


DO 


LD A, M(X+) 


X NOW POINTS TO C-EXP. 


282 F6A4 


D2 


LD A, M(X-) 




283 F6A5 


9C00 


IFEQ A, 




284 F6A7 


B4F974 


JMPL UNDFL 




285 F6AA 


90FE 


IFGT A, OFE 




286 F6AC 


B4F980 


JMPL OVRFL 




287 F6AF 


F2 


LD A, W(X-) 




288 F6B0 


F2 


LD A, W(X-) 


X NOW POINTS TO C-LO. 


289 F6B1 


B5F9C8 + 


JSR FPAK 


PACK C. 


290 F6B4 


3FC4 


POP SP 


SET UP SP FOR RETURN. 


291 F6B6 


3FCC 


POP B 




292 F6B8 


3FCE 


POP X 




293 F6BA 


3C 


RET 




294 




; 




295 




.END 
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35 






.FORM 'FMULT.MAC' 


36 






.INCLD FMULT.MAC 


1 






.TITLE FMULT 


2 
3 

4 






.LOCAL 




SUBROUTINE TO MULTIPLY TWO SP FLOATING POINT NUMBERS. 


5 
6 
7 






C = F1*F2 




Fl IS 


STORED IN THE IEEE FORMAT IN REGS K AND A. 


8 




THE HIGH WORD OF Fl WILL BE REFERRED AS Fl-Rl AND IS IN K. 


9 




THE LOW WORD OF Fl WILL BE REFERRED TO AS Fl-RO AND IS IN A. 


10 








11 




F2 IS 


STORED IN THE IEEE FORMAT ON THE STACK. IF SP IS THE 


12 




STACK 


POINTER ON ENTRY, THEN 


13 




THE HIGH WORD OF F2 , REFERRED TO AS F2-R1 IS AT SP - 4 AND 


14 




THE LOW WORD OF F2, REFERRED TO AS F2-R0 IS AT SP - 6. 


15 








16 




C IS 


RETURNED IN THE IEEE FORMAT IN REGS K AND A. 


17 




REGS. 


X AND B ARE PRESERVED. 


18 








19 


I 


MULT: 




20 


F6BB AFCE 




PUSH X ; SAVE X ON ENTRY. 


21 


F6BD AFCC 




PUSH B ; SAVE B ON ENTRY. 


22 




SAVE 


4DDRESS OF F2-R0 IN IMP1. 


23 


F6BF ACC4CE 




LD X, SP 


24 


F6C2 86FFF6CEF8 




ADD X, 0FFF6 ; SUBTRACT 10. 


25 


F6C7 ACCEOO 




LD TMP1, X ; SAVE IN TMP1. 


26 




CHECK 


AND SEE IF Fl IS A NAN. 


27 


F6CA B5F973 + 




JSR FNACHK 


28 


F6CD 07 




IF C 


29 


F6CE B4F93E 




JMPL FNAN ; Fl IS A NAN. 


30 




CHECK 


AND SEE IF F2 IS A NAN. 


31 


F6D1 ACCACC 




LD B, K 


32 


F6D4 ACC8CE 




LD X, A 


33 


F6D7 A20200A8 




LD A, W(TMPl+2) 


34 


F6DB ACC8CA 




LD K, A 


35 


F6DE AECE 




X A, X 


36 


F6E0 B5F95D + 




JSR FNACHK 


37 


F6E3 07 




IF C 


38 


F6E4 B4F928 




JMPL FNAN ; F2 IS NAN. 


39 




CHECK 


AND SEE IF F2 IS ZERO. 


40 


F6E7 B5F967 + 




JSR FZCHK 


41 


F6EA 07 




IF C 


42 


F6EB 94DC 




JMP JCZERO ; F2 IS ZERO. 


43 




CHECK 


AN SEE IF Fl IS ZERO. 


44 


F6ED ACCCCA 




LD K, B ; RESTORE Fl-Rl FROM B. 


45 


F6FD B5F95E + 




JSR FZCHK 


46 


F6F3 07 




IF C 


47 


F6F4 94D3 




JMP JCZERO ; Fl IS ZERO. 


48 




GET HERE MEANS NORMAL MULTIPLICATION. 


49 




UNPACK Fl AND F2. 
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50 F6F6 


ACC4CE 


51 F6F9 


8210C4F8 


52 F6FD 


AFCE 


53 F6FF 


B5F95F 


54 F702 


ACOOCC 


55 F705 


ACCEOO 


56 F708 


EO 


57 F709 


40 


58 F70A 


AECA 


59 F70C 


E4 


60 F70D 


AECA 


61 F70F 


B5F94F 


62 




63 F712 


F2 


64 F713 


ACOOCC 


65 F716 


E2 


66 F717 


40 


67 




68 F718 


F4 


69 F719 


C7 


70 F71A 


ACC8CA 


71 F71D 


E4 


72 F71E 


B9FF00 


73 F721 


C7 


74 F722 


96CAF8 


75 F725 


F6 


76 F726 


E2 


77 F727 


40 


78 F728 


99FF 


79 F72A 


C7 


80 F72B 


8FFB 


81 F72D 


B8C080 


82 F730 


07 


83 F731 


46 


84 




85 F732 


E7 


86 F733 


07 


87 F734 


B4F8E7 


88 F737 


C7 


89 




90 




91 F738 


E7 


92 F739 


07 


93 F73A 


B4F8F2 


94 F73D 


96C817 


95 F740 


96C808 


96 F743 


F3 


97 




98 




99 




100 F744 


F2 



LD X, SP 
ADD SP, 010 
PUSH X 
JSR FUNPAK 
LD B, TMP1 
LD TMP1, X 
LOS A, W(B+) 
NOP 

X A, K 
LD A, W(B) 
X A, K 
JSR FUNPAK 

SET X TO POINT TO 
LD A, W(X-) 
LD B, TMP1 
LDS A, W(B-) 
NOP 

COMPUTE C-EXP AND 
W(X) 



X POINTS TO Fl-LO. 

MOVE SP PAST LOCAL STORAGE. 

SAVE SP ON STACK FOR QUICK RETURN. 

UNPACK Fl. 

B NOW POINTS TO F2-R0 . 

TMP1 NOW POINTS TO F2-L0. 

LOAD F2-R0 INTO A. 



; LOAD F2-R1 INTO K. 
; UNPAK F2. 
F2-SIGN AND B TO POINT TO Fl-SIGN. 



LD A 
SHR A 
LD K, 
LD A, 
AND A 
SHR A 
ADD A, K 
ST A, W(X) 



A 

W(B) 
OFFOO 



LOS A, W(B-) 

NOP 

AND A, OFF 

SHR A 

XOR A, W(X) 

ADD A, 0C080 

IF C 

JP $EXCH2 
; CHECK TO SEE IF EXP IS ZERO 

SHL A 

IF C 

JMPL UNDFL 

SHR A 

; CHECK FOR EXPONENT OVERFLOW 
$EXCH2 : 

SHL A 

IF C 

JMPL OVRFL 

IF A. 7 

SET A.O 

X A, W(X-) 



SIGN AND STORE IN F2-EXP AND F2-SIGN. 
A IS (EEEEEEEE-F2) . (SSSSSSSS-F2) 
SHR SINCE SUM OF EXPS CAN BE 9 BITS. 
K IS (DEEEEEEEE-F2) . (SSSSSSS-F2) 
A IS (EEEEEEEE-F1) . 
MASK OUT SIGN BITS. 
A IS (DEEEEEEEE-F1) . 
A IS (EEEEEEEEE-C) . 
STORE IN F2-SIGN. 
A IS (EEEEEEEE-F1) . 



(SSSSSSSS-F1) 

,(0000000) 

(SSSSSSS-F2) 

(SSSSSSSS-F1) 



MASK OUT EXP BITS. 

A IS (000000000SSSSSSS-F1) 

A IS (EEEEEEEEESSSSSSS-C) 

REMOVE EXCESS BIAS OF 127-DEC FROM EXP. 



IF CARRY, THEN NO UNDERFLOW NOW 
IF NOT, UNDERFLOW FOR SURE. 



UNDERFLOW, SO JUMP. 

RESTORE BIT SHIFTED OUT (0). 



IF C IS 1, 

THEN OVERFLOW FOR SURE. 



RESTORE LAST BIT OF SIGN. 

STORE C-EXP. C-SIGN IN F2-EXP.F2-SIGN. 



MULTIPLY THE MANTISSA. 
FIRST COMPUTE F1-HI*F2-HI. 
LD A, W(X-) 
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101 


F745 


ACCEOO 


102 


F748 


FE 


103 


F749 


A6FFFAC4AB 


104 


F74E 


AECE 


105 


F750 


A6FFFCC4AB 


106 






107 


F755 


ACOOCE 


108 


F758 


FO 


109 


F759 


ACCEOO 


110 


F75C 


FE 


111 


F75D 


AECE 


112 


F75F 


A6FFFAC4F8 


113 


F764 


A6FFFAC4AB 


114 


F769 


07 


115 


F76A 


A6FFFCC4A9 


116 






117 


F76F 


E2 


118 


F770 


40 


119 


F771 


ACOOCE 


120 


F774 


F4 


121 


F775 


FE 


122 


F776 


AECE 


123 


F778 


A6FFFAC4F8 


124 


F77D 


A6FFFAC4AB 


125 


F782 


A6FFFCC4AB 


126 


F787 


07 


127 


F788 


04 


128 






129 






130 


F789 


ACOOCE 


131 


F78C 


BD7FFF 


132 


F78F 


4D 


133 






134 


F790 


E7 


135 


F791 


F3 


136 


F792 


A6FFFAC4AB 


137 


F797 


E7 


138 


F798 


Fl 


139 


F799 


07 


140 


F79A 


8F08 


141 


F79C 


51 


142 






143 






144 


F79D 


F3 


145 


F79E 


A6FFFAC4A8 


146 


F7A3 


Fl 


147 


F7A4 


FO 


148 


F7A5 


F4 


149 


F7A6 


B80100 


150 


F7A9 


07 


151 


F7AA 


B4F882 



TMP1 NOW POINTS TO F2-L0 . 

STORE LOW WORD OF PRODUCT ON STACK. 

STORE HIGH WORD OF PRODUCT ON STACK. 

TMP1 NOW POINTS TO F2-HI . 

ADD LOW WORD OF LAST PROD. TO HIGH WORD. 



LD TMP1, X ; 

MULT A, W(B) 

ST A, W(SP-6) ; 

X A, X 

ST A, W(SP-4) ; 
; NOW COMPUTE F1-HI*F2-L0 

LD X, TMP1 

LD A, W(X+) 

LD TMP1, X ; 

MULT A, W(B) 

X A, X 

ADD A, W(SP-6) ; 

ST A, W(SP-6) 

IF C 

INC W(SP-4) ; IF CARRY, INCREASE HIGH WORD BY 1. 
; FINALLY COMPUTE F1-L0*F2-HI. 

LDS A, W(B-) ; ADJUST B TO POINT TO Fl-LO. 

NOP 

LD X, TMP1 

LD A, W(X) 

MULT A, W(B) 

X A, X 

ADD A, W(SP-6) 

ST A, W(SP-6) 

LD A, W(SP-4) 

IF C 

INC A 



ADD LOW WORD ACCUMULATED SO FAR. 

A CONTAINS HIGH WORD OF PRODUCT. 
IF CARRY ON LAST LOW WORD ADD, 
THEN INCREASE HIGH WORD. 



; MANTISSA MULTIPLICATION DONE. NOW CHECK FOR NORMALIZATION. 
LD X, TMP1 

IFGT A, 07FFF ; IS MSB OF PRODUCT 1 ? 
JP $EXINC ; YES, INCREASE MANTISSA. 

; NEED TO SHIFT MANTISSA LEFT BY 1 BIT. 
SHL A 
X A, W(X-) 
LD A, W(SP-6) 
SHL A 
X A, W(X+) 
IF C 

SET W(X) .0 
JP {ROUND 



DID SHIFT OF LOW WORD PUSH OUT A 1 ? 

YES SO SET LSB OF HIGH WORD. 

GO TO ROUNDING CODE. 
JEXINC : 
; NEED TO INCREASE EXPONENT BY 1. REMEMBER X POINTS TO F2-HI. 



X A, W(X-) 
LD A, W(SP-6) 
X A, W(X+) 
LD A, W(X+) 
LD A, W(X) 
ADD A, 0100 
IF C 
JMPL OVRFL 



A CONTAINS HIGH WORD, X POINTS TO F2-L0. 
GET LOW WORD. 
STORE LOW WORD. 

GET C-EXP.C-SIGN 
INCREASE C-EXP. 

EXPONENT OVERFLOW. 
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152 F7AD 


F3 




X A, W(X-) ; 


NO OVERFLOW, SO SAVE 


153 




; ROUNDING CODE. 




154 




$ROUND : 






155 F7AE 


B5F8ED 




JSRL SROUND 




156 




; FINAL 


CHECK OF EXPONENT 


. 


157 F7B1 


DO 




LD A, M(X+) ; 


X NOW POINTS TO C-EX 


158 F7B2 


D2 




LD A, M(X-) 




159 F7B3 


9C00 




IFEQ A, 




160 F7B5 


B4F866 




JMPL UNDFL 




161 F7B8 


9DFE 




IFGT A, OFE 




162 F7BA 


B4F872 




JMPL OVRFL 




163 F7BD 


F2 




LD A, W(X-) 




164 F7BE 


F2 




LD A, W(X-) ; 


X NOW POINTS TO C-LO 


165 F7BF 


B5F8BA 


+ 


JSR FPAK ; 


PACK C. 


166 F7C2 


3FC4 




POP SP ; 


SET UP SP FOR RETURN 


167 F7C4 


3FCC 




POP B 




168 F7C6 


3FCE 




POP X 




169 F7C8 


3C 




RET 




170 




; EXCEPTION HANDLING. 




171 




; C IS 


ZERO B'COS ONE OF 


Fl OR F2 IS ZERO. 


172 




$CZERO: 






173 F7C9 


00 




CLR A 




174 F7CA 


ACC8CA 




LD K, A 




175 F7CD 


3FCC 




POP B 




176 F7CF 


3FCE 




POP X 




177 F7D1 


3C 




RET 




178 




; 






179 






.END 





.C-SIGN. 
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37 




.FORM 'FDI V.MAC' 




38 




.INCLD FDI V.MAC 




1 




.TITLE FDIV 




2 
3 

4 




.LOCAL 




SUBROUTINE TO DIVIDE TWO SP FLOATING POINT NUMBERS. 




5 
6 
7 


C 


= F1/F2 




Fl IS 


STORED IN THE IEEE FORMAT IN REGS K AND A. 




8 


THE HIGH WORD OF Fl WILL BE REFERRED AS Fl-Rl AND IS IN K. 


9 


THE LOW WORD OF Fl WILL BE REFERRED TO AS Fl-RO AND IS 


IN A. 


10 








11 


F2 IS 


STORED IN THE IEEE FORMAT ON THE STACK. IF SP IS 


THE 


12 


STACK 


POINTER ON ENTRY, THEN 




13 


THE HIGH WORD OF F2 , REFERRED TO AS F2-R1 IS AT SP - 4 


AND 


14 


THE LOW WORD OF F2, REFERRED TO AS F2-R0 IS AT SP - 6. 




15 








16 


C IS 


RETURNED IN THE IEEE FORMAT IN REGS K AND A. 




17 








18 I 


DIV: 






19 F7D2 AFCE 




PUSH X 




20 F7D4 AFCC 




PUSH B 




21 


SAVE 


4DDRESS OF F2-R0 IN TMP1. 




22 F7D6 ACC4CE 




LD X, SP 




23 F7D9 86FFF6CEF8 




ADD X, 0FFF6 ; SUBTRACT 10. 




24 F7DE ACCEOO 




LD TMP1, X ; AND SAVE IN TMP1. 




25 


CHECK 


AND SEE IF Fl IS A NAN. 




26 F7E1 85F85C + 




JSR FNACHK 




27 F7E4 07 




IF C 




28 F7E5 B4F827 




JMPL FNAN ; Fl IS A NAN. 




29 


CHECK 


AND SEE IF F2 IS A NAN. 




30 F7E8 ACCACC 




LD B, K 




31 F7EB ACC8CE 




LD X, A 




32 F7EE A20200A8 




LD A, W(TMPl+2) 




33 F7F2 ACC8CA 




LD K, A 




34 F7F5 AECE 




X A. X 




35 F7F7 B5F846 + 




JSR FNACHK 




36 F7FA 07 




IF C 




37 F7FB B4F811 




JMPL FNAN ; F2 IS NAN. 




38 


CHECK 


AND SEE IF F2 IS ZERO. 




39 F7FE B5F850 + 




JSR FZCHK 




40 F801 07 




IF C 




41 F802 B4F7FB 




JMPL DIVBYO ; F2 IS ZERO. 




42 


CHECK 


AND SEE IF Fl IS ZERO. 




43 F805 ACCCCA 




LD K, B ; RESTORE Fl-Rl FROM B. 




44 F808 B5F846 + 




JSR FZCHK 




45 F80B 07 




IF C 




46 F80C 94F1 




JMP $CZERO ; Fl IS ZERO. 




47 


GET HERE MEANS NORMAL DIVISION. 




48 


UNPACK Fl AND F2. 




49 F80E ACC4CE 




LD X, SP ; X POINTS TO Fl-LO. 
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50 F811 


8210C4F8 


51 F815 


AFCE 


52 F817 


B5F847 


53 F81A 


ACOOCC 


54 F81D 


ACCEOO 


55 F820 


EO 


56 F821 


40 


57 F822 


AECA 


58 F824 


E4 


59 F825 


AECA 


60 F827 


B5F837 


61 




62 




63 




64 F82A 


F2 


65 F82B 


F2 


66 F82C 


ACOOCC 


67 F82F 


E2 


68 F830 


40 


69 F831 


E2 


70 F832 


40 


71 F833 


ACC8CA 


72 F836 


F4 


73 F837 


FD 


74 F838 


51 


75 




76 




77 F839 


EO 


78 F83A 


40 


79 F83B 


AECA 


80 F83D 


B80100 


81 F840 


E3 


82 F841 


40 


83 F842 


E4 


84 F843 


C7 


85 F844 


E3 


86 F845 


40 


87 F846 


E4 


88 F847 


D7 


89 F848 


El 


90 F849 


40 


91 




92 




93 




94 F84A 


FO 


95 F84B 


EO 


96 F84C 


40 


97 F84D 


F4 


98 F84E 


B9FF00 


99 F851 


C7 


100 F852 


ACC8CA 



ADD SP, 010 

PUSH X 

JSR FUNPAK 

LD B, TMP1 

LD TMP1, X 

LDS A, W(B+) 

NOP 

X A, K 

LD A, W(B) 

X A, K 

JSR FUNPAK 



MOVE SP PAST LOCAL STORAGE. 

SAVE SP ON STACK FOR QUICK RETURN. 

UNPACK Fl. 

B NOW POINTS TO F2-R0 

TMP1 NOW POINTS TO F2-L0. 

LOAD F2-R0 INTO A. 



LOAD F2-R1 INTO K. 
UNPAK F2. 



ENSURE THAT Fl-HI IS LESS THAN F2-HI . 



LD A, 


W(X-) ; 


LD A, 


W(X-) ; 


LD B, 


TMP1 ; 


LDS A 


, W(B-) ; 


NOP 




LDS A 


, W(B-) ; 


NOP 




LD K, 


a ; 


LD A, 


W(X) ; 


IFGT 


1, W(B) ; 


JP JFEXSN ; 


LOS A 


, W(B+) ; 


NOP 




X A, 


£ ; 


ADD A 


, 0100 ; 


XS A, 


W(B-) ; 


NOP 




LD A, 


W(B) ; 


SHR A 




XS A, 


W(B-) ; 


NOP 




LD A, 


W(B) ; 


RRC A 




XS A, 


W(B+) ; 


NOP 





X POINTS TO F2-EXP.F2-SIGN. 

X POINTS TO F2-HI. 

B POINTS TO F2-L0. 

B POINTS TO F1-EXP.F1-SIGN. 

LOAD F1-EXP.F1-SIGN. 

B POINTS TO Fl-HI. 

SAVE F1-EXP.F1-SIGN IN K. 

LOAD F2-HI. 

IS F2-HI > Fl-HI ? 

YES, SO ALL IS WELL. 

GET HERE MEANS NEED TO SHR Fl, 

AND INCREASE ITS EXPONENT. 

GET Fl-HI. 

B POINTS TO F1-EXP.F1-SIGN. 

SWAP F1-EXP.F1-SIGN AND Fl-HI. 

INCREASE Fl-EXP BY 1. 

STORE BACK IN Fl-EXP. Fl-SIGN. 

B POINTS TO Fl-HI. 

LOAD Fl-HI. 

STORE BACK IN Fl-HI. 
B POINTS TO Fl-LO. 
LOAD Fl-LO. 

PUT II BACK IN Fl-LO. 
B POINTS TO Fl-HI. 



$FEXSN : 
; DETERMINE C-EXP AND C-SIGN. 



LD A, 
LDS A, 
NOP 
LD A, 
AND A, 
SHR A 
LD K, 



W(X+) 
W(B+) 

W(X) 
OFFOO 



X POINTS TO F2-EXP.F2-SIGN. 
B POINTS TO Fl-EXP. Fl-SIGN. 

LOAD F2-EXP.F2-SIGN. 

MASK OUT THE SIGN. 

ALLOW 9 BITS FOR EXP CALCULATIONS. 

SAVE IT IN K. 
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101 


F855 


E4 


102 


F856 


B9FF00 


103 


F859 


C7 


104 


F85A 


02 


105 


F85B 


96CAEB 


106 






107 






108 


F85E 


B7000000 


109 


F862 


E7 


110 


F863 


07 


111 


F864 


B700FF00 


112 


F868 


D7 


113 


F869 


B83F00 


114 


F86C 


E7 


115 


F86D 


06 


116 


F86E 


49 


117 


F86F 


960010 


118 


F872 


B4F7A9 


119 


F875 


B4F7B7 


120 






121 






122 






123 






124 


F878 


BCFFOO 


125 


F87B 


B4F7B1 


126 


F87E 


9C00 


127 


F880 


B4F79B 


128 


F883 


ABOO 


129 


F885 


F4 


130 


F886 


99FF 


131 


F888 


FB 


132 


F889 


99FF 


133 


F88B 


9600FA 


134 


F88E 


F3 


135 






136 


F88F 


F2 


137 


F890 


E2 


138 


F891 


40 


139 






140 






141 






142 


F892 


FO 


143 


F893 


ACCEOO 


144 


F896 


FE 


145 






146 






147 


F897 


AECC 


148 


F899 


AEOO 


149 


F89B 


AECC 


150 






151 


F890 


EF 



LD A, W(B) 
AND A, OFFOO 
SHR A 
SET C 
SUBC A, K 



LD TMP1, 

SHL A 

IF C 

LD TMP1, 

RRC A 

ADD A, 03F00 

SHL A 

IFN C 

JP $FSIGN 

IF TMP1.0 

JMPL UNDFL 

JMPL OVRFL 



OFF 



LOAD F1-EXP.F1-SIGN. 
MASK OUT SIGN. 



SUBTRACT THE EXPONENTS. 

NOTE THAT NOW THE MS 9 BITS 

OF A CONTAIN A 2'S COMP. INTEGER. 



SAVE SIGN OF NUMBER IN TMP1. 

RESTORE IEEE BIAS. 

MAKE EXPONENT 8 BITS. 

NO CARRY ? 

THEN ALL IS WELL. 

WAS EXP NEGATIVE BEFORE ? 

YES, SO UNDERFLOW. 

OTHERWISE OVERFLOW. 



$FSIGN: 

; C-EXP HAS BEEN COMPUTED. NOW FIND C-SIGN. 
; BUT FIRST TAKE CARE OF SPECIAL OVER/UNDERFLOW CASES. 
IFEQ A, OFFOO 



SAVE C-EXP. 00000000 IN TMP1. 

LOAD F2-EXP.F2-SIGN. 

MASK OUT F2-EXP. 

A NOW HAS Fl-EXP. C-SIGN. 

MASK OUT Fl-EXP. 

BRING IN C-EXP. 

STORE IN F2-EXP.F2-SIGN. 

X POINTS TO F2-HI. 

X POINTS TO F2-L0. 

B POINTS TO Fl-HI. 



JMPL 


OVRFL 


IFEQ 


A, 


JMPL 


UNDFL 


ST A, 


TMP1 


LD A, 


W(X) 


AND A 


, OFF 


XOR A 


, W(B) 


AND A 


, OFF 


OR A, 


TMP1 


X A, 


W(X-) 


LD A, 


W(X-) 


LDS A 


, W(B- 


NOP 





NOW DO THE MANTISSA DIVISION. 



LD A, W(X+) 
LD TMP1, X 
MULT A, W(B) 



X A, B 
X A, TMP1 
X A, B 

.BYTE OEF 



LOAD F2-L0. X POINTS TO F2-HI . 
SAVE ADDRESS OF F2-HI IN TMP1. 
COMPUTE F2-L0*F1-HI. 
X CONTAINS MS WORD AND A IS LS WORD. 



A POINTS TO Fl-HI, 
A POINTS TO F2-HI, 
A CONTAINS LS WORD. 



B CONTAINS LS WORD. 
TMP1 POINTS TO Fl-HI. 
B POINTS TO F2-HI. 



DIVD A, W(B) 



KLUDGED 
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153 F89E 


ACC8CA 


LD K, A 


SAVE QUOTIENT IN K. 


154 F8A1 


A8CC 


LD A, B 


A POINTS TO F2-HI. 


155 F8A3 


AEOO 


X A, TMP1 


A POINTS TO Fl-HI, TMP1 POINTS TO F2-HI . 


156 F8A5 


AECC 


X A, B 


B POINTS TO Fl-HI. 


157 F8A7 


E2 


LDS A, W(B-) 


B POINTS TO Fl-LO. 


158 F8A8 


40 


NOP 




159 F8A9 


EO 


LOS A, W(B+) 


LOAD Fl-LO. 


160 F8AA 


40 


NOP 


B POINTS TO Fl-HI. 


161 F8AB 


02 


SET C 




162 F8AC 


96CAEB 


SUBC A, K 


SUBTRACT QUOTIENT SAVED IN K. 


163 F8AF 


ACC8CE 


LD X, A 


AND SAVE IN X. 


164 F8B2 


E4 


LD A, W(B) 


LOAD Fl-HI. 


165 F8B3 


06 


IFN C 


IF C WAS NOT SET IN THE LAST SUBTRACT, 


166 F8B4 


05 


DEC A 


ADJUST THE BORROW. 


167 F8B5 


AECE 


X A, X 




168 F8B7 


ACOOCC 


LD B, TMP1 


B POINTS TO F2-HI. 


169 




; 




170 F8BA 


EF 


.BYTE OEF 


DIVD A, W(B) - KLUDGED AGAIN ! 


171 






QUOTIENT IN A, REM IN X. 


172 




; 




173 F8BB 


ABOO 


ST A, TMP1 


SAVE QUOTIENT IN TMP1. 


174 F8BD 


00 


CLR A 


ZERO A. 


175 




; 




176 F8BE 


EF 


.BYTE OEF 


DIVD A, W(B) - KLUDGED YET AGAIN ! 


177 




; 




178 F8BF 


AEOO 


X A, TMP1 


SWAP OLD AND NEW QUOTIENTS. 


179 




; 




180 




; CHECK FOR NORMALIZATION. CAN BE OFF BY AT MOST 1 BIT. 


181 F8C1 


E7 


SHL A 




182 F8C2 


07 


IF C 




183 F8C3 


56 


JP JNMED 


IT IS NORMALIZED. 


184 






GET HERE MEANS NEED TO SHIFT LEFT ONCE. 


185 F8C4 


AEOO 


X A, TMP1 


SWAP HI AND LO WORDS. 


186 F8C6 


E7 


SHL A 




187 F8C7 


AEOO 


X A, IMP1 


HI WORD IS IN A, LO WORD IN TMP1. 


188 F8C9 


07 


IF C 


WAS 1 SHIFTED OUT OF LO WORD? 


189 F8CA 


96C808 


SET A.O 


YES, THEN SET LSB OF HI WORD. 


190 F8CD 


ABC A 


ST A, K 


SAVE HI WORD IN K. 


191 F8CF 


El 


XS A, W(B+) 




192 F8D0 


40 


NOP 


B POINTS TO F2-EXP.F2-SIGN. 


193 F8D1 


E4 


LD A, W(B) 


LOAD F2-EXP.F2-SIGN. 


194 F8D2 


B8FF00 


ADD A, OFFOO 


SUBTRACT 1 FROM EXPONENT. 


195 F8D5 


E3 


XS A, W(B-) 


STORE BACK IN F2-EXP.F2-SIGN. 


196 F8D6 


40 


NOP 


B POINTS TO F2-HI. 


197 F8D7 


A8CA 


LD A, K 


HI WORD TO A. 


198 F8D9 


E7 


SHL A 




199 




$NMED : 




200 F8DA 


D7 


RRC A 


RESTORE BIT OUT. 


201 F8DB 


E3 


XS A, W(B-) 


SAVE HI-WORD IN F2-HI . 


202 F8DC 


40 


NOP 


B POINTS TO F2-L0. 
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203 F8DD 


A800 




LD A, TMP1 






204 F8DF 


El 




XS A, W(B+) 




; SAVE C-LO. 


205 F8E0 


40 




NOP 




; B POINTS TO F2-HI . 


206 F8E1 


ACCCCE 




LD X, B 




; MOVE ADDRESS OF F2-HI TO X 


207 




; 








208 




; ROUNDING CODE. 






209 F8E4 


B5F7B7 




JSRL SROUND 






210 




; FINAL 


CHECK OF EXPONENT. 


211 F8E7 


DO 




LD A, M(X+) 




; X NOW POINTS TO C-EXP. 


212 F8E8 


D2 




LD A, M(X-) 






213 F8E9 


9C00 




IFEQ A, 






214 F8EB 


B4F730 




JMPL UNDFL 






215 F8EE 


9DFE 




IFGT A, OFE 






216 F8F0 


B4F73C 




JMPL OVRFL 






217 F8F3 


F2 




LD A, W(X-) 






218 F8F4 


F2 




LD A, W(X-) 




; X NOW POINTS TO C-LO. 


219 F8F5 


B5F784 


+ 


JSR FPAK 




; PACK C. 


220 F8F8 


3FC4 




POP SP 




; SET UP SP FOR RETURN. 


221 F8FA 


3FCC 




POP B 






222 F8FC 


3FCE 




POP X 






223 F8FE 


3C 




RET 






224 




; C IS 


ZERO B'COS Fl 


IS 


ZERO. 


225 




$CZERO : 








226 F8FF 


00 




CLR A 






227 F900 


ACC8CA 




LD K, A 






228 F903 


3FCC 




POP B 






229 F905 


3FCE 




POP X 






230 F907 


3C 




RET 






231 




; 








232 






.END 
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.FORM 'FSINX.MAC 


40 
1 
2 






.INCLD FSINX.MAC 






.TITLE SINX 


3 






.LOCAL 




4 








A VERY DIRTY APPROXIMATION TO SIN(X). 


5 
6 
7 








X SHOULD BE IN RADIANS. 










ON INPUT X SHOULD BE IN IEEE FLP FORMAT IN REGS. K AND A. 


8 

9 

10 








ON RETURN SIN(X) IS IN 


IEEE FLP FORMAT IN REGS. K AND A. 






SINX: 




11 


F908 


AFCE 


PUSH X 


SAVE X. 


12 


F90A 


AFC8 


PUSH A 




13 


F90C 


AFC A 


PUSH K 


X TO THE STACK. 


14 


F90E 


3653 


JSRL FMULT 


COMPUTE X»2. 


15 


F910 


AFC8 


PUSH A 




16 


F912 


AFC A 


PUSH K 


X A 2 TO THE STACK. 


17 


F914 


B6F994A8 


LD A, W($A5L0) 




18 


F918 


A4F996CAAB 


LD K, W($A5HI) 


LOAD A5. 


19 


F91D 


3662 


JSRL FMULT 


COMPUTE A5*X A 2. 


20 


F91F 


AFC8 


PUSH A 




21 


F921 


AFC A 


PUSH K 




22 


F923 


B6F998A8 


LD A, W($A4L0) 




23 


F927 


A4F99ACAAB 


LD K, W($A4HI) 


LOAD A4. 


24 


F92C 


B5FBE6 


JSRL FSUB 


COMPUTE A4-A5*X A 2. 


25 


F92F 


3FCE 


POP X 




26 


F931 


3FCE 


POP X 




27 


F933 


3678 


JSRL FMULT 


COMPUTE 


28 








X A 2(A4 - A5*X A 2) . 


29 


F935 


AFC8 


PUSH A 




30 


F937 


AFC A 


PUSH K 




31 


F939 


B6F99CA8 


LD A, W($A3L0) 




32 


F93D 


A4F99ECAAB 


LD K, W($A3HI) 


LOAD A3. 


33 


F942 


B5FBD0 


JSRL FSUB 


COMPUTE 


34 








A3 - X A 2(A4 - A5*X A 2) . 


35 


F945 


3FCE 


POP X 




36 


F947 


3FCE 


POP X 




37 


F949 


368E 


JSRL FMULT 


COMPUTE 


38 








X A 2(A3 - X*2(A4 - A5*X A 2) ) . 


39 


F94B 


AFC8 


PUSH A 




40 


F94D 


AFC A 


PUSH K 




41 


F94F 


B6F9A0A8 


LD A, W($A2L0) 




42 


F953 


A4F9A2CAAB 


LD K, W($A2HI) 


LOAD A2. 


43 


F958 


B5FBBA 


JSRL FSUB 


COMPUTE 


44 








A2 - X A 2(A3 - X*2(A4 - A5*X*2) ) . 


45 


F95B 


3FCE 


POP X 




46 


F95D 


3FCE 


POP X 




47 


F95F 


36A4 


JSRL FMULT 


COMPUTE 


48 








X A 2(A2 - X A 2(A3 - X A 2(A4 - A5*X A 2) ) ) . 


49 


F961 


AFC8 


PUSH A 
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50 


F963 


AFCA 




PUSH K 




51 


F965 


B6F9A4AB 




LD A, W($A1L0) 




52 


F969 


A4F9A6CAAB 




LD K, W($A1HI) 


LOAD Al. 


53 


F96E 


B5FBA4 




JSRL FSUB 


COMPUTE 


54 










Al - X A 2(A2 - X*2(A3 - X A 2(A4 - A5*X A 2) ) ) . 


55 


F971 


3FCE 




POP X 




56 


F973 


3FCE 




POP X 




57 


F975 


36BA 


- 


JSRL FMULT 


COMPUTE 


58 










X*2(A1 - X*2(A2 - X A 2(A3 - X*2(A4 - A5*X A 2) ) ) ) . 


59 


F977 


AFC8 




PUSH A 




60 


F979 


AFCA 




PUSH K 




61 


F97B 


B13F80 




LD K, 03F80 




62 


F97E 


00 




CLR A 


LOAD 1.0 INTO K-A. 


63 


F97F 


B5FB93 




JSRL FSUB 


COMPUTE 


64 










1 - ALL THE JUNK ABOVE. 


65 


F982 


3FCE 




POP X 




66 


F984 


3FCE 




POP X 




67 


F986 


3FCE 




POP X 




68 


F988 


3FCE 




POP X 


NOW X IS AT THE TOP OF STACK. 


69 


F98A 


36CF 


- 


JSRL FMULT 


COMPUTE 


70 








; X(l - X A 2(A1 - X*2(A2 - X A 2(A3 - X*2(A4 - A5*X*2) ) ) ) ) . 


71 


F98C 


3FCE 




POP X 




72 


F98E 


3FCE 




POP X 




73 


F990 


3FCE 




POP X 




74 


F992 


3C 




RET 




75 






; 






76 


F993 


40 




.EVEN 




77 






; 






78 


F994 


2B32 


$A5L0 


.WORD 0322B 




79 


F996 


D732 


$A5HI 


.WORD 032D7 




80 


F998 


1DEF 


JA4L0 


.WORD 0EF1D 




81 


F99A 


3836 


JA4HI 


.WORD 03638 




82 


F99C 


010D 


$A3L0 


.WORD 00D01 




83 


F99E 


5039 


$A3HI 


.WORD 03950 




84 


F9A0 


8988 


$A2L0 


.WORD 08889 




85 


F9A2 


083C 


$A2HI 


.WORD 03C08 




86 


F9A4 


ADAA 


$A1L0 


.WORD OAAAD 




87 


F9A6 


2A3E 


$A1HI 


.WORD 03E2A 




88 












89 






; A DIRTY APPROXIMATION TO COS(X) USING SIN(X). 


90 












91 






COSX: 






92 


F9A8 


AFCE 




PUSH X 




93 


F9AA 


ACC8CE 




LD X, A 




94 


F9AD 


B6F9C8A8 




LD A, W($PI2L0) 




95 


F9B1 


AFC8 




PUSH A 




96 


F9B3 


B6F9CAA8 




LD A, W($PI2HI) 




97 


F9B7 


AFC8 




PUSH A 




98 


F9B9 


A8CE 




LD A, X 




99 


F9BB 


B5FB77 




JSRL FADD 


COMPUTE X + PI/2. 


100 


F9BE 


3FCE 




POP X 
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101 


F9C0 


3FCE 


102 


F9C2 


34BA 


103 


F9C4 


3FCE 


104 


F9C6 


3C 


105 






106 


F9C7 


40 


107 


F9C8 


DBOF 


108 


F9CA 


C93F 


109 






110 






111 






112 






113 


F9CC 


AFCE 


114 


F9CE 


AFCC 


115 


F9D0 


AFC8 


116 


F9D2 


AFC A 


117 


F9D4 


342 C 


118 


F9D6 


ACC8CE 


119 


F9D9 


ACCACC 


120 


F9DC 


3FCA 


121 


F9DE 


3FC8 


122 


F9E0 


AFCE 


123 


F9E2 


AFCC 


124 


F9E4 


34DC 


125 


F9E6 


3614 


126 


F9E8 


3FCC 


127 


F9EA 


3FCC 


128 


F9EC 


3FCC 


129 


F9EE 


3FCE 


130 


F9F0 


3C 


131 






132 






41 






42 






43 


FFFE 


OOFO 



COMPUTE SIN(X+PI/2) . 



POP X 
JSRL SINX 
POP X 
RET 

.EVEN 
$PI2L0: .WORD OOFDB 
$PI2HI: .WORD 03FC9 



A DIRTY APPROXIMATION TO TAN(X) USING SINX AND COSX. 

TANX: 

PUSH X 
PUSH B 
PUSH A 
PUSH K 
JSR COSX 
LD X, A 
LD B, K 
POP K 
POP A 
PUSH X 
PUSH B 
JSR SINX 
JSR FDIV 
POP B 
POP B 
POP B 
POP X 
RET 

.END 



COMPUTE COS(X) 



COMPUTE SIN(X) . 
COMPUTE TAN(X) 



SIN(X)/COS(X) . 



.END LISTER 
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A 


00C8 W 


AT OF 


F13A * 


B 


OOCC W 


BFMUL 


F0C7 


cosx 


F9A8 


DIVBYO 


FOOO 


FADD 


F535 


FDIV 


F7D2 


FMULT 


F6BB 


FNACHK 


F040 


FNAN 


FOOF 


FPAK 


F07C 


FPERWD 


0002 W 


FPTRAP 


F09D 


FSUB 


F515 


FTOA 


F311 * 


FUNPAK 


F061 


FZCHK 


F051 


ISIOK 


F105 


K 


OOCA W 


LISTER 


FOOO 


MUL10 


F118 


OVRFL 


F02F 


PC 


00C6 W 


SINX 


F908 


SP 


00C4 W 


SROUND 


F09E 


TANX 


F9CC * 


TMP1 


0000 w 


UNDFL 


F01E 


X 


OOCE W 


$A10EX 


F1F4 


$A1HI 


F9A6 


$A1L0 


F9A4 


$A2HI 


F9A2 


$A2L0 


F9A0 


$A3HI 


F99E 


$A3L0 


F99C 


$A4HI 


F99A 


}A4L0 


F998 


$A5HI 


F996 


$A5L0 


F994 


$ACCF 


F1A9 


$ACCM 


F177 


$ADDEX 


F1FF 


$ADDMN 


F612 


$ADEM 


F3EC 


$ADJEX 


F685 


$ANORM 


F64A 


$ANOTO 


F05C 


$CHK0T 


F113 


pHNGS 


F361 


$CSIGN 


F349 


$CZERO 


F7C9 


$CZERO 


F8FF 


$DIV10 


F272 


pIVIO 


F370 


$DOLUP 


F481 


$DTHI 


F270 


$DTHI 


F392 


$DTLO 


F26E 


pTLO 


F390 


JESAVE 


F20F 


$ESIGN 


F1C7 


pXACC 


F1C9 


$EXCH2 


F738 


$EXCHR 


F1BB 


$EXCLP 


F1D4 


pXCOL 


F1BA 


PXCPT 


F2CF 


$EXIN2 


F0B7 


$EXINC 


F79D 


((EXIT 


F0C6 


$F1CHK F56E 


$F1SIN 


F643 


$F2GTR 


F5EB 


$TEXSN 


F84A 


$FRCOL 


F18B 


$FSIGN 


F878 


pOBAK 


F41F 


$G00N 


F42B 


$HIUP 


FOAF 


JINCOL 


F15C 


$INCRV 


F40C 


$INDUN 


F445 


$INTFY 


F43C 


$ISNAN 


F04C 


$ISNED 


F298 


$ISNED 


F3BE 


$ISNXT 


F17D 


JJAMDN 


F29E 


$JAMDN 


F3C4 


$JAMIT 


F280 


$JAMIT 


F3A6 


JJAMLP 


F287 


$JAMLP 


F3AD 


$L00P1 


F147 


$L00P1 


F5FB 


$L00P2 


F5D7 


$ML4 


F3EA 


$ML0G2 


F33F 


$MSIGN 


F154 


$MTHI 


F26C 


$MTHI 


F396 


$MTLO 


F26A 


$MTL0 


F394 


$MUL10 


F398 


pAGAS 


F2C3 


pAN 


F4CF 


pANLP 


F4DC 


$NEG10 


F20B 


pLOOP 


F653 


pMED 


F8DA 


P0RM1 


F230 


$N0RM2 


F235 


pRDUN 


F247 


pRLUP 


F237 


pTZER 


F585 


$0V1 


F666 


$0VR1 


F29B 


$0VR1 


F3C1 


$PI2HI 


F9CA 


$PI2L0 


F9C8 


PEMV9 


F354 


$RNDUP 


F0A7 


POUND 


F6A0 


$ROUND 


F7AE 


$TRADD 


F66C 


$VDOWN 


F418 


$VUP 


F416 


$ZERLP 


F505 


$ZERO 


F4E9 


$ZR0F1 


F607 


$ZR0F2 


F5E2 











52 



NATIONAL SEMICONDUCTOR CORPORATION PAGE: 50 

HPC CROSS ASSEMBLER, REV:C, 30 JUL 86 

SINX 

MACRO TABLE 

NO WARNING LINES 

NO ERROR LINES 

2547 ROM BYTES USED 

SOURCE CHECKSUM = A31F 
OBJECT CHECKSUM = 2AC3 

INPUT FILE C:LISTER.MAC 
LISTING FILE C:LISTER.PRN 
OBJECT FILE C:LISTER.LM 
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LIFE SUPPORT POLICY 



NATIONAL'S PRODUCTS ARE NOT AUTHORIZED FOR USE AS CRITICAL COMPONENTS IN LIFE SUPPORT 
DEVICES OR SYSTEMS WITHOUT THE EXPRESS WRITTEN APPROVAL OF THE PRESIDENT OF NATIONAL 
SEMICONDUCTOR CORPORATION. As used herein: 



1. Life support devices or systems are devices or 2. A critical component is any component of a 



life 



systems which, (a) are intended for surgical implant 
into the body, or (b) support or sustain life, and whose 
failure to perform, when properly used in accordance 
with instructions for use provided in the labeling, can 
be reasonably expected to result in a significant injury 
to the user. 



support device or system whose failure to perform can 
be reasonably expected to cause the failure of the life 
support device or system, or to affect its safety or 
effectiveness. 
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National Semiconductor 
Corporation 

1111 West Bardin Road 
Arlington, TX 76017 
Tel: 1(800) 272-9959 
Fax: 1(800) 737-7018 



National Semiconductor 
Europe 

Fax: (+49) 0-180-530 85 86 

Email: cnjwge@tevm2.nsc.com 

Deutsch Tel: ( + 49) 0-180-530 85 85 

English Tel: (+49) 0-180-532 78 32 

Frangais Tel: ( + 49) 0-180-532 93 58 

Italiano Tel: (+49) 0-180-534 16 80 



National Semiconductor 
Hong Kong Ltd. 

13th Floor, Straight Block, 
Ocean Centre, 5 Canton Rd. 
Tsimshatsui, Kowloon 
Hong Kong 
Tel: (852)2737-1600 
Fax: (852) 2736-9960 



National Semiconductor 
Japan Ltd. 

Tel: 81-043-299-2309 
Fax: 81-043-299-2408 



National does not assume any responsibility for use of any circuitry described, no circuit patent licenses are implied and National reserves the right at any time without notice to change said circuitry and specifications. 



